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ABSTRACT 


This publication contains the preliminary agenda and summaries for the Third Annual JPL 
Airborne Geoscience Workshop, held at the Jet Propulsion Laboratory, Pasadena, California, on 
June 1—5, 1992. This main workshop is divided into three smaller workshops as follows: 

• The Airborne Visible/Infrared Imaging Spectrometer (AYIRIS) workshop, on 
June 1 and 2. The summaries for this workshop appear in Volume 1. 

. The Thermal Infrared Multispectral Scanner (TIMS) workshop, on June 3. The 
summaries for this workshop appear in Volume 2. 

• The Airborne Synthetic Aperture Radar (AIRSAR) workshop, on June 4 and 5. 
The summaries for this workshop appear in Volume 3. 
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12:00 noon Lunch 

1:00 p.m. Mapping Vegetation Types With the Multiple Spectral Feature Mapping 
Algorithm in Both Emission and Absorption 

Roger N. Clark, Gregg A. Swayze, Christopher Koch, and Cathy Ager 
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1:30 p.m. 
2:00 p.m. 

2:30 p.m. 

3:00 p.m. 
3:30 p.m. 

4:00 p.m. 

4:30 p.m. 

5:00 p.m. 
5:30 p.m. 
5:45 p.m. 


Agenda (continued) 

Third annual jpl Airborne Geoscience workshop: 

Airborne visible/infrared imaging Spectrometer 

(AVTRIS) 

Multiple Dataset Water-Quality Analyses in the Vicinity of an Ocean Wastewater 
Plume 

Michael Hamilton, Curtiss O. Davis, W. Joseph Rhea, and 
Jeannette van den Bosch 

MAC Europe 91: Evaluation of AVIRIS, GER Imaging Spectrometry Data for 
the Land Application Testsite Oberpfaffenhofen 

F. Lehmann, R. Richter, H. Rothfuss, K. Werner, P. Hausknecht, A. Muller, and 
P. Strobl 

Using Endmembers in AVIRIS Images To Estimate Changes in Vegetative 
Biomass 

Milton O. Smith, John B. Adams, Susan L. Ustin, and Dar A. Roberts 
Break 

Atmospheric Correction of AVIRIS Data in Ocean Waters 
Gregory Terrie and Robert Arnone 

The 1991 AVIRIS/POLDER Experiment in Camargue, France 
F . Baret, C, Lepneur , S. Jacquemoud, V. Carr&re, X.F. Gu, M. Steven, 

V. Vanderbilt, J F . Hanocq, S. Ustin, G. Rondeaux, C. Daughtry, L. Biehl, 

R. Pettigrew, D. Modro, H. Horoyan, T. Sarto, C. Despontin, and 
H. Razafindraibe 

AVIRIS Sensor and Ground Data System: Status and Plans 
Thomas Chrien and Earl Hansen 

Wrap up. 

End of AVIRIS Workshop. 

Shuttle bus departs JPL for the Pasadena Ritz-Carlton Hotel. 
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AGENDA 


THIRD ANNUAL JPL AIRBORNE GEOSCIENCE WORKSHOP: 

THERMAL INFRARED MULTISPECTRAL SCANNER 
(TIMS) 

June 3, 1992 
Von Karman Auditorium 
Jet Propulsion Laboratory 
Pasadena, California 91109 

WEDNESDAY, JUNE 3, 1992 

7:15 a.m. Shuttle bus departs Pasadena Ritz-Carlton Hotel for JPL. 

7:30 a.m. Registration and continental breakfast at JPL. 

8:00 a.m. Welcome 

8:30 a.m. TIMS Performance Evaluation Summary 

Bruce Spiering, G. Meeks , J. Anderson , S. Juggi, o.nd S> Kuo 

A Quantitative Analysis of TIMS Data Obtained on the Leaijet 23 at Various 

Altitudes 

S. Jaggi 

Analysis of TIMS Performance Subjected to Simulated Wind Blast 
S. Jaggi and S. Kuo 

9:30 a.m. Sensitivity of Blackbody Reference Panels to Wind Blast 

Gordon Hoover 

10:00 a.m. Break 

10 30 a.m. Preliminary Analysis of TIMS Performance on the ER-2 
S.J. Hook, VJ. Realmuto, and R.E. Alley 

1 1:00 a.m. Comparison of Preliminary Results From Airborne ASTER Simulator (A AS) 

With TIMS Data . ^ _ 

Yoshiaki Kannari, Franklin Mills, Hiroshi Watanabe, Teruya Ezaka, 

Tatsuhiko Narita, and Sheng-Huei Chang 

1 1 :30 a.m. Simulation of ASTER Data Using AVIRIS Images 

Michael Abrams 

12:00 noon Lunch 

1 :00 p.m. Application of Split Window Technique to TIMS Data 

Tsuneo Matsunaga, Shuichi Rokugawa, and Yoshinori Ishii 

1:30 p.m. Atmospheric Corrections for TIMS Estimated Emittance 
TA. Warner and D.W. Levandowski 

2:00 p.m. An Algorithm for the Estimation of Bounds on the Emissivity and Temperatures 
From Thermal Multispectral Airborne Remotely Sensed Data 
S. Jaggi, D. Quattrochi, andR. Baskin 
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AGENDA (CONTINUED) 


2:30 p.m. 

3:00 p.m. 
3:30 p.m. 

4:00 p.m. 
4:30 p.m. 

5:00 p.m. 

5:30 p.m. 
5:45 p.m. 


THIRD ANNUAL JPL AIRBORNE GEOSCIENCE WORKSHOP: 

Thermal Infrared multispectral Scanner 
(TIMS) 

Multi-Resolution Processing for Fractal Analysis of Airborne Remotely Sensed 
Data 

S. Jaggi, D. Quattrochi, and N. Lam 
Break 

Preliminary Analysis of Thermal-Infrared Multispectral Scanner Data of the Iron 
Hill, Colorado Carbonatite-Alkalic Rock Complex 
Lawrence C. Rowan, Kenneth Watson, and Susanne H. Miller 

The Use of TIMS for Mapping Different Pahoehoe Surfaces: Mauna Iki, Kilauea 
Scott K. Rowland 

Ejecta Patterns of Meteor Crater, Arizona Derived From the Linear Un-Mixing of 
TIMS Data and Laboratory Thermal Emission Spectra 
Michael S. Ramsey and Philip R. Christensen 

The Use of TIMS Data To Estimate the SO 2 Concentrations of Volcanic Plumes: 
A Case Study at Mount Etna, Sicily 
Vincent J. Realmuto 

End of TIMS Workshop. 

Shuttle bus departs JPL for the Pasadena Ritz-Carlton Hotel. 
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AGENDA 


THIRD ANNUAL JPL AIRBORNE GEOSCIENCE WORKSHOP: 
AIRBORNE SYNTHETIC APERTURE RADAR 
(AIRSAR) 

June 4 and 5, 1992 
Von Karman Auditorium 
Jet Propulsion Laboratory 
Pasadena, California 91109 


THURSDAY, 
7:15 a.m. 
7:30 a.m. 
8:00 a.m. 
8:30 a.m. 

9:00 a.m. 

9:30 a.m. 


10:00 a.m. 
10:30 a.m. 

11:00 a.m. 

11:30 a.m. 

12:00 noon 
1:00 p.m. 

1:30 p.m. 

2:00 p.m. 

2:30 p.m. 


JUNE 4, 1992 

Shuttle bus departs Pasadena Ritz-Carlton Hotel for JPL. 


Registration and continental breakfast at JPU. 


Welcome 


The NASA/JPL Three-Frequency AIRSAR System 
J. van Zyl, R. Carande, Y. Low, T. Miller, and K. Wheeler 


A Snow Wetness Retrieval Algorithm for SAR 
Jiancheng Shi and Jeff Dozier 


Comparison of JPL- AIRSAR and DLR E-SAR Images from the MAC Europe ’91 
Campaign Over Testsite Oberpfaffenhofen: Frequency and Polarization 
Dependent Backscatter Variations From Agricultural Fields 
C. Schmullius and J. Nithack 


Break 


Monitoring Environmental State of Alaskan Forests With AIRSAR 

Kyle C. McDonald, JoBea Way, Eric Rignot, Cindy Williams, Les Viereck, and 

Phylis Adams 


Comparison of Modeled Backscatter With JAR Data at P-Band 
Yong Wang, Frank W. Davis, and John M. Melack 


SAR Backscatter From Coniferous Forest paps 
John L. Day and Frank W. Davis 


Lunch 


Panel Discussion on Future Emphasis of the AIRSAR System 
J. van Zyl, Moderator 


Retrieval of Pine Forest Biomass Using JPL AIRSAR Data 
A. Beaudoin, T. Le Loan, F. Zagolski, C. C Hsu. H.C. Han, andJA. Kong 

Characterization of Wetland, Forest, and Agricultural Ecosystems in Belize With 

Airborne Radar (AIRSAR) . 

Kevin O. Pope, Jose Maria Rey-Benayas, and Jack t . Pans 


Strategies for Detection of Floodplain Inundation With Multi-Frequency 
Polarimetric SAR 

Laura L. Hess and John M. Melack _ 
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Agenda (continued) 


Third Annual jpl airborne Geoscience workshop: 
Airborne Synthetic Aperture Radar 
(AIRSAR) 


3:00 p.m. 
3:30 p.m. 

4:00 p.m. 

4:30 p.m. 
5:00 p.m. 
5:30 p.m. 
5:45 p.m. 
6:30 p.m. 
9:00 p.m. 


Break 

Supervised Fully Polarimetric Classification of the Black Forest Test Site: From 

MAESTRO 1 to MAC Europe 

G. De Grandi, C. Lavalle, H. De Groof, and A. Sieber 

Relating Multifrequency Radar Backscattering to Forest Biomass: Modeling and 
AIRSAR Measurement 
Guoqing Sun and K. Jon Ranson 

Poster Previews 

Poster Previews 

End of session. 

Shuttle bus departs JPL for the Pasadena Ritz-Carlton Hotel. 

Reception and poster sessions at the Pasadena Ritz-Carlton Hotel. 

Close of reception and poster sessions. 
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AGENDA 


THIRD ANNUAL JPL AIRBORNE GEOSCIENCE WORKSHOP: 

poster Session 

Thursday, June 4, 1992 
6:30 to 9:00 p.m. 

Pasadena Ritz-Carlton Hotel 

Processing of AIRSAR Polarimetric Data for Soil Moisture Estimation Over 
Mahantango Watershed Area 

K. S. Rao 

Evaluation of Polarimetric S AR Parameters for Soil Moisture Retrieval 
Jiancheng Shi, Jakob J. van Zyl, and Edwin T. Engman 

Interaction Types and Their Like-Polarization Phase-Angle Difference Signatures 
Jack F. Paris 

Application of Modified VICAR/IBIS GIS to Analysis of July 1991 Flevoland 
AIRSAR Data 

L. Norikane, B. Broek, and A. Freeman 

Radar Analysis and Visualization Environment (RAVEN): Software for 

Polarimetric Radar Analysis 

K.S. Kierein-Young, A.B. Lefkoff, and FA. Kruse 

Measuring Ocean Coherence Time With Dual-Baseline Interferometry 
Richard E. Carande 

A Bibliography of Global Change, Airborne Science, 1985-1991 
Edwin J. Sheffner and James G. Lawless 

ATTIRE (Analytical Tools for Thermal Infrared Engineering) — A Thermal 
Sensor Simulation Package 
5. Jaggi 

Kilauea Data Set Complied for Distribution on Compact Disc 

Lori Glaze, George Karas, Sonia Chernobieff, Elsa Abbott, and Earnie Paylor 

The JPL Spectral Library 0.4 to 2.5 Micrometers 
Simon J. Hook, Cindy /. Grove, and Earnest D. Paylor II 

Lossless Compression of AVIRIS Data: Comparison of Methods and Instrument 
Constraints 

R.E. Roger, J.F. Arnold, M.C. Cavenor, andJ.A. Richards 

Simulation of AVHRR-K Band Ratios With AVIRIS 
Melanie A. Wetzel and Ronald M. Welch 

Soil Conservation Applications With C-Band SAR 
B. Brisco, RJ. Brown, J. Naunheimer, and D. Bedard 
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Agenda 


Third Annual JPL Airborne Geoscience workshop: 

Thursday Poster Session (continued) 

14. Comparison of Edges Detected at Different Polarisations in MAESTRO Data 
Ronald G. Caves, Peter J. Harley, and Shaun Quegan 

15. Identification of Erosion Hazards in a Mediterranean Environment 
M. Altherr, J. Hill, and W. Mehl 
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AGENDA 


THIRD ANNUAL JPL AIRBORNE GEOSCIENCE WORKSHOP: 
AIRBORNE SYNTHETIC APERTURE RADAR 
(AIRSAR) 

June 4 and 5, 1992 
Von Kantian Auditorium 
Jet Propulsion Laboratory 
Pasadena, California 91109 


FRIDAY, JUNE 5, 1992 


7:15 a.m. 
7:30 a.m. 
8:00 a.m. 

8:30 a.m. 
9:00 a.m. 


Shuttle bus departs Pasadena Ritz-Carlton Hotel for JPL. 


Registration and continental breakfast at JPL. 


Oceanic Features Detected by SAR in the Mediterranean Sea During the MAC 
Europe ’91 Campaign 
Werner Alpers 

SAR Observations in the Gulf of Mexico 
David Sheres 


Investigation of AIRSAR Signatures of the Gulf Stream ■ v u r 

G.R. Valenzuela, J.S. Lee, D.L. Schuler, G.O. Marmorino, F. Askari, K. Hoppel, 
J.A.C. Kaiser, and W.C. Keller 


9:30 a.m. 

10:00 a.m. 
10:30 a.m. 


11:00 a.m. 


11:30 a.m. 


Mapping of Sea Bottom Topography 

C.J. Calkoen, GJ. Wensink, and G.H.F.M. Hesselmans 

Break 

Sea Bottom Topography Imaging With SAR 
M.WA. van derKooij, GJ. Wensink, andJ. Vogelzang 

AIRSAR Surveys of Upper-Ocean Fronts Off California and Hawaii 
P. Flament 

Preliminary Results of Polarization Signatures for Glacial Moraines in the Mono 
Basin, Eastern Sierra Nevada 

Richard R. Forster, Andrew N. Fox, and Bryan hacks 


12:00 noon Lunch 

1 :00 p.m. Detecting Surface Roughness Effects on the Atmospheric Boundary Layer Via 
AIRSAR Data: A Field Experiment in Death Valley, California 
Dan G. Blumberg and Ronald Greeley 

1 :30 p.m. Extraction of Quantitative Surface Characteristics From AIRSAR Data for Death 
Valley, California 
K.S. Kierein-Young and FA. Kruse 


2:00 p.m. The TOPSAR Interferometric Radar Topographic Mapping Instrument 
Howard A. Zebker, Spren N. Madsen, Jan Martin, Giovanni Alberti, 
Sergio Vetrella, and Alessandro Cucci 
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Agenda (continued) 


2:30 p.m. 

3:00 p.m. 
3:30 p.m. 

4:00 p.m. 

4:30 p.m. 

5:00 p.m. 
5:30 p.m. 
5:45 p.m. 


Third Annual JPL Airborne Geoscience workshop: 

Airborne Synthetic Aperture Radar 
(AIRSAR) 

Evaluation of the TOPSAR Performance by Using Passive and Active Calibrators 
G. Alberti, A. Moccia, S. Ponte, and S. Vetrella 

Break 

Fitting a Three-Component Scattering Model to Polarimetric SAR Data 
A. Freeman and S. Durden 

Application of Symmetry Properties to Polarimetric Remote Sensing With JPL 
AIRSAR Data 

S.V. Nghiem, S.H. Yueh, R. Kwok, and F.K. Li 

External Calibration of Polarimetric Radar Images Using Distributed Targets 
Simon H. Yueh, S.V. Nghiem, and R. Kwok 

Wrap up. 

End of AIRSAR Workshop. 

Shuttle bus departs JPL for the Pasadena Ritz-Carlton Hotel. 
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A SNOW WETNESS RETRIEVAL ALGORITHM FOR SAR 

Jiancheng Shi 

Center for Remote Sensing and Environmental Optics 
University of California, Santa Barbara, CA 93106 
JefT Dozier 

NASA Godddard Space Flight Center, Code 900, 

Greenbelt, MD 20771 


The objectives of this study are (1) to evaluate the backscattering signals re- 
sponse to snow wetness and (2) to develop an algorithm for snow wetness measure- 
ment using C-band polarimetric SAR. In hydrological investigations, modeling and 
*3 forecasting of snowmelt runoff requires information about snowpack properties and 

their spatial variability. In particular, timely measurement of snow parameters is 
needed for operational hydrology. The liquid water content of snowpack is one of 
the important parameters. Active microwave sensors are highly sensitive to liquid 
water in the snowpack because of the large dielectric contrast between ice and water 
in the microwave spectrum. They are not affected by weather and have a spatial 
resolution compatible with the topographic variation in alpine regions. However, a 
quantitative algorithm for retrieval snow wetness has not been developed. 

The data, from a polarimetric SAR overflight by the NASA/JPL DC-8 over 
the glaciers in the Oztal Alps, Austria, were analyzed. In order to utilize the po- 
larimetric SAR data for quantitatively measuring snow wetness purpose, we must 
consider the accuracy of the radar measurements. For radiometric calibration, and 
as geometric reference points, trihedral Conner reflectors with leg lengths of 0.65, 
1.25 and 1.8 meters were deployed on the glaciers. Using topographic information, 

~ SAR data can be calibrated (van Zyl and Shi, 1992) so that it contains the more 

accurate information of the polarization signatures response to snow physical pa- 
rameters. With the generated incidence angle and radar antenna angle map, the 
topographic effect on SAR image data is minimized so that a direct comparison in 
the model simulations and inversion snow wetness become possible. 

Backscattering measurements by SAR from wet snow covered terrain are af- 
fected by two sets of parameters (Shi and Dozier, 1992a): (1) sensor parameters 
which include the frequency, polarization, and viewing geometry, and (2) snow- 
pack parameters which include snow density, liquid water content, particle sizes 
and shapes of ice and water, type of the correlation function and its parameters of 
surface roughness. 

In addition to snowpack physical parameters, the surface roughness and in- 
cidence angle have a great impact on the relationship between the backscattering 
coefficients and snow wetness. A negative correlation between the backscattering 
coefficients and snow wetness was found by the truck-mounted Microwave Active 
Spectrometer at frequency 8.6 GHz for incidence angles greater than 5° (Stiles and 
Ulaby, 1980). However, positive correlation was found in this study when compar- 
ing AIRSAR measurements of low and high snow wetness samples. The observed 
backscattering coefficients for VV and IIIT polarizations at identical incidence angles 
were about 5 dB higher from high liquid water content snow samples than from low 
liquid water content snow samples. A 3 dB difference was found in IIV polarization. 

Both measurements and model predictions show that the scattering mecha- 
nisms of wet snow are characterized as (1) for low liquid water content snow, surface 
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scattering is dominant at small incidence angles and volume scattering dominates 
at large incidence angles. (2) for high liquid water content snow with a relatively 
rough surface, surface scattering dominates at all incidence angles. 

The actual relationship between the co-polarization signals and snow wetness 
is controlled by the scattering mechanisms* When the surface is smooth, volume 
scattering is the dominant scattering source. As snow wetness increases, both the 
volume scattering albedo and the transmission coefficients greatly decrease. This re- 
sults in a negative correlation between the backscattering signals and snow wetness. 
When the surface is not smooth, increasing snow wetness results in greatly increased 
surface scattering interaction and surface scattering becomes the dominant scatter- 
ing process. Therefore, a positive correlation between the backscattering signals 
and snow wetness will be observed. 

To develop an algorithm of measuring snow wetness, we need to minimize or 
combine these factors since the polarimetric SAR only provides a limited number 
of independent observations. The task is to select appropriate measurements for 
relating snow wetness to the backscattering measurements. 

Based on the identified scattering mechanics of wet snow-covered terrain from 
the model predictions and measurements of the polarimetric properties (Shi and 
Dozier, 1992a), we construct the inversion model with two components: 


<T? P = *7 


+ < 


pp 


( 1 ) 


where pp indicates polarization. a t is the total backscattering coefficient. <j s is the 
surface backscattering from the air-snow interface and a v is the volume backscat- 
tering from the snowpack. 

Since the relations of the ratios derived by the first-order surface and volume 
scattering model in co-polarization signals, two ratios of the combined VV and IT II 
polarization signals can be represented as 


T hh 
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where Dr and Dt are the surface and volume backscattering ratios of VV and IIH 
polarizations (Shi and Dozier, 1992b), which are only a function of incidence angle 
and the permittivity of wet snow. The Chh is Hie surface backscattering contribution 
in the ITII polarization of the total backscattering signal. From the above two 
measurements, the two unknowns, e r and Chh , ca » be solved. The algorithm derived 
above does not require any information about the surface roughness and the volume 
scattering albedo. It only involves the calculation of snowpack permittivity. 

To test the algorithm for measuring snow wetness over a large area, a snow map 
was first obtained by performing the supervised Bayes classification (Shi et al., 1991) 
and non-snow-covered area was masked. Secondly, the backscattering coefficients 
of VV and IIII polarizations for a given pixel were determined by the median value 
within a 5 x 5 window to reduce the efTect of image speckle. The algorithm was 
applied after the measurements of a hh /(a vv + cr hh ) and a vv /(cr vv cr hh ) were tested 
by the physical conditions (Shi and Dozier, 1992b). 

Figure shows an image of the inversion-derived snow wetness. The image 
brightness is proportional to the snow wetness and ranges from 1 to 8 percent 
by volume. The black region is non-snow- covered area. At the time of the radar 
survey the snow cover was wet at all elevation zones. The liquid water content of 
the top snow layer was in the order of 5 to 6 percent by volume at the elevation 
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of 3,000 meters, decreasing to 2 or 3 percent at the highest elevations. The snow 
densities and depths ranged from 492 to 580 kg m -3 and from 45 to 114 cm re- 
spectively. Using co-registered DEM data, we compared the inversion- derived snow 
wetness within the elevation zones. The result shows that the algorithm performs 
well and that the absolute error, generally, is less than 2 percent. 

This study shows recent results of our efforts to develop and verify an algo- 
rithm for snow wetness retrieval from a polarimetric SAR. Our algorithm is based 
on the first-order scattering model with consideration of both surface and volume 
scattering. It operates at C-band and requires only rough information about the ice 
volume fraction in snowpack. By evaluating the relationship between the backscat- 
tering coefficients of both surface and volume scattering for the co-polarization sig- 
nals, the estimate of snow wetness can be derived from the ratios of the combined 
co- polarization signals. The inversion algorithm performs well using AIRSAR data 
and should prove useful for routine and large-area snow wetness measurements. 



Figure Inversion-derived snow wetness. The image brightness is proportional to 
the snow wetness and ranges from 1 to 8 percent by volume. The black region is 
nonsnow-covered area. 


References 
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COMPARISON OF JPL-AIRSAR AND DLR E-SAR IMAGES 
FROM THE MAC EUROPE '91 CAMPAIGN OVER TESTSITE OBER- 
PFAFFENHOFEN: FREQUENCY AND POLARIZATION DEPENDENT 
BACKSCATTER VARIATIONS FROM AGRICULTURAL FIELDS 

C. Schmullius, J. Nithack _ T m 

N » 4 

Institute of Radio-Frequency Technology 
German Aerospace Establishment (DLR) 

D-8031 Oberpfaffenhofen, Germany 

1. INTRODUCTION 

On July 12, the MAC Europe '91 (Multi-Sensor Airborne 
Campaign) took place over test site Oberpfaffenhofen. The 
DLR Institute of Radio-Frequency Technology participated 
with its C-W, X-W, and X-HH Experimental SAR (E-SAR). The 
high resolution E-SAR images with a pixel size between 1 
and 2 m and the polarimetric AIRSAR images were analyzed. 
Using both sensors in combination is a unique opportunity 
to evaluate SAR images in a frequency range from p- to 
X-band and to investigate polarimetric information. 

2. EXPERIMENT AND DATA DESCRIPTION 

The JPL AIRSAR was flown over testsite Oberpfaffenhofen on 
July 12, 1991 at noon. The DLR E-SAR was flown on different 
daytimes between July 8 and 12. The E-SAR data were proces- 
sed at the DLR facility with the Motion Compensation SAR 
Processor [ Moreira 1990 ] . The pixel resolution for the 
8-look images (slant range presentation) is 1.8 x 1.2 m. 
The antenna pattern was corrected using a multiplicative 
polynomial correction function. The image digital numbers 
correspond to 16 bit amplitude values. 

Ground truth consisted in mapping of landuse and the pheno- 
logical situation, and intensive soil sampling on large, 
flat agricultural fields belonging to three different soil 
types: a loessy soil, drained turfy material, and a gra- 

vel rich glacial terrace. Biomass and leaf area index were 
measured on several corn fields. 

3. IMAGE INTERPRETATION 

Digital image processing was performed on an EASI/PACE sy- 
stem (PCI Inc.) using a DECstation 5000/200. Since the E— 
SAR samples one frequency-polarization combination at a 
time, the image channels of the test site had to be geome- 
trically rectified to each other. Figure 1 shows a color 
composite of three E-SAR images illustrating different 
backscattering mechanisms for C- and X-band but, furthermo- 
re, also for X-HH and -W polarization. These differences 
are due to varying states of crop maturity. The wetness, 
mainly in the crop's ears, causes absorption losses of the 
W-polarized signals. 
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This finding is in agreement with measurements of the ab- 
sorption loss factor of wheat stalks [Allen and Ulaby 1984] 
and negative correlations between biomass and C-W back- 
scatter for corn fields and potato fields f see figure 2. 
For vertical polarization the loss factor increases with 
plant volumetric moisture (ie. the backscattered power 
decreases) , and for horizontal polarization the loss 
factor decreases with wetness. Hence, the HH-polarized 
returns penetrate the crop's ear layer somewhat and back- 
scattering occurs on lower parts of the plants where the 
different state of maturity does not cause, at this time of 
the year (July 8), differences in water content. 

Figure 3 shows an unsupervised Migrating-Means (or 
K-Cluster) approach to classify an agricultural area close 
to the DLR. Only the E-SAR C-W and X-W channels were 
used. Image processing techniques such as 
post-classification filtering or a combination with image 
segmentation algorithms [Rignot et al. 1991 ] can improve 
the classification results. 

4. CONCLUSION AND OUTLOOK 

The W-polarized C- and X-band scattering is dependent on 
the different degrees of maturity in the barley fields. 
HH-polarized X-band signals have a larger penetration depth 
and the backscattering intensity is not dependent on the 
water content of the uppermost vegetation layer. SAR images 
of one week later do not show these backscattering variati- 
ons since the crops had been ripening in that time period. 

jPL polarimetric AIRSAR data from this test site is being 
analysed with special regard to the polarization-dependent 
backscatter phenomena, since it could be shown that a com- 
bination of a different frequency plus a different polari- 
zation gives best statistical separability. Future cam- 
paigns such as CLEOPATRA 1992 will include additional 
ground truth measurements that are regarded to be of impor- 
tance and have been missing in the MAC Europe '91 Campaign, 
such as plant geometry. 

5. REFERENCES 

C. Allen and F. Ulaby, "Modeling the Polarization Dependen- 
ce of the Attenuation in Vegetation Canopies," in Procee - 
d ings of the International Geoscience and Remote Sensing 
Symposium , Strasbourg, France, 27.-30.8. 1984. 

J. Moreira, "Motion Compensation SAR-Processing Facility at 
DLR, " in Proceedings of the 10th EARSeL Symposium , Toulou- 
se, 5.-8. June 1990, pp. 279-288. 

E. Rignot, R. Chellappa, R. Dubois, R. Kwok, and J. van 
Zyl, "Unsupervised Segmentation of Polarimetric SAR Data 
using the Covariance Matrix," in Proceedings of the Inter - 
national Geoscience and Remote Sensing Symposium , Espoo, 
3.-6. June 1991, pp. 1813-1816. 
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Fig. 1 

Color composite of a three channel image 
(red: C-W, green: X-W, blue: X-HH, MAC 
Europe '91). 


correlation <t. y fresh cn □ n ass 



with the fresh biomass of com and potatoes 
(C-band image, Aug. 1989). 

Fig. 3: 

Non-supervised classification of a two fre- 
quency image (X-band, C-band both 
W-polarisation, MAC Europe '91). 
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During March 1988 and May 1991, the JPL airborne synthetic aperture radar, 
AIRSAR, collected sets of multi-temporal imagery of the Bonanza Creek Experi- 
mental Forest near Fairbanks, Alaska. These data sets consist of series of multi- 
polarized images collected at P-, L- and C-bands each over a period of a few days. 
The AIRSAR campaigns were complemented with extensive ground measurements 
that included observations of both static canopy characteristics such as forest ar- 
chitecture as well as properties that vary on short term time scales such as canopy 
dielectric conditions. Observations exist for several stands of deciduous and conif- 
erous species including white spruce ( Picea glauca), black spruce (Picea martana) 
and balsam poplar ( Populus balsamifera). 

Although the duration of each campaign was fairly short, significant changes in 
environmental conditions caused notable variations in the physiological state of the 
canopies. During the 1988 campaign, environmental conditions ranged from unsea- 
sonably warm to more normal subfreezing temperatures. This permitted AIRSAR 
observations of frozen and thawed canopy states. During May 1991, ice jams that 
occurred along the river caused many stands to flood while the subsequent clear- 
ing of the river then allowed the waters to recede, leaving a snow covered ground 
surface. This allowed observations of several stands during both flooded and non- 
flooded conditions. Furthermore, the local weather varied from clear sunny days 
to heavy overcast days with some occurrence of rain. Measurements of leaf water 
potential indicated that this caused significant variations in canopy water status, 
allowing SAR observations of water stressed and unstressed trees. 

Mean backscatter from several stands is examined for the various canopy physi- 
ological states. The changes in canopy backscatter that occur as a function of envi- 
ronmental and physiological state are analyzed. Preliminary results of a backscatter 
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signature modeling analysis are presented. The implications of using SAR to mon- 
itor canopy phenological state are addressed. 


This work was carried out at the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract to the National Aeronautics and Space Administration, and at the 
Institute of Northern Forestry, University of Alaska. 
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1. Introduction 

In recent years several analytical models have been developed to predict 
microwave scattering by trees and forest canopies [1-6]. These models contribute to the 
understanding of radar backscattcr over forested regions to the extent that they capture 
the basic interactions between microwave radiation and tree canopies, understories, and 
ground layers as functions of incidence angle, wavelength, and polarization. The Santa 
Barbara microwave model backscatter model for woodland (i.e. with discontinuous tree 
canopies) combines a single-tree backscattcr model and a gap probability model [5-6]. 
Comparison of model predictions with synthetic aperture radar (SAR) data at L-band (A, 
= 0.235 m) is promising [7], but much work is still needed to test the validity of model 
predictions at other wavelengths. Here we test the validity of the model predictions at 
P-band (A = 0.68 m) for woodland stands at our Mt. Shasta test site. 


2. Study area 

Extensive ground data have been collected to support our SIR-B (Shuttle Imaging 
Radar) and SIR-C/X-SAR studies. The site is located to the southeast of Mt. Shasta, 
California (41°18’ N, 122°05’ W), and spans elevations from 1160 to 1220 m. The 
forest stands under investigation are natural, are on level ground, and are dominated by 
pondcrosa pine (Pinus ponderosa) or pine mixed with white fir (Abies concolor). Under 
story vegetation is sparse and consists primarily of perennial grasses and Jorbes. The 
litter layer can reach a depth of about 0.1 m. Three ponderosa pine stands ( Sp2 , St2 , 
and "Sll 1") in the site were chosen for this study because adequate ground data for these 
stands are available to provide inputs to the model. The three stands have different «and 
densities and trunk diameters at breast height (dbh). "Sp2" has trees of 25.0 m h and 
0.77 m dbh, with stand density of 23 trees per hectare. "Stll" has trees of 27.0 m h and 
0.46 m dbh, with stand density of 228 trees per hectare. "St2" has trees of 25 m h and 
0 42 m dbh, with stand density of 303 trees per hectare. The soil of the stands is derived 
from recent alluvial deposits of volcanic ash, and the ground surface is smooth at P- 
band. Other ground data used as model inputs can be found in [8]. 

3. SAR data 

SAR data were acquired by the Jet Propulsion Laboratory (JPL) airborne SAR 
overflights on 6 September 1989 after a dry summer. The SAR data were processed and 
calibrated by JPL. The estimated calibration uncertainty of the backscatter is ±1.0 dB 
[3] The estimated calibration uncertainty of VV-HH phase difference is ±10°, and the 
estimated calibration uncertainty of HH and VV correlation coefficient is ±0.1 (personal 
communication. Freeman 1992). We received the standard 4-look compressed data with 
pixel spacing of 12.1 m (azimuth) and 6.7 m (slant range). For the three stands, the 
means of the backscatter (HH, HV, and VV polarizations), the VV-HH phase difference, 
and the HH and VV correlation coefficient were extracted from the SAR data. 
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4. Comparison of model results with SAR data at P-band 

Model simulations for stands "St2", "Stl 1”, and "Sp2" were carried out at the same 
incidence angle (0J as the SAR data. We have 8 SAR data takes covering stand "St2" 
(9 0 ranges from 22° to 54°), 5 data takes covering stand "Stl 1" (Q 0 from 25° to 45°), and 
5 data takes covering stand , ’Sp2" (0* from 27° to 47°). The HH, HV, and VV back- 
scatter, and the HH and VV correlation coefficient are shown in Figure 1. The model 
makes good predictions of the HH backscatter for stands "St2" and "Stll" (Figure la), 
the VV backscatter for all three stands (Figure lb), and the HV backscatter for stand 
"Stl 1" (Figure lc). The model overestimates the HH backscatter for stand "Sp2" (Figure 
la). The model underestimates the HV backscatter for stand "Sp2" and stand "St2" (Fig- 
ure lc). We attribute the underestimate to the surface model, in which we model only 
the co-polarizcd surface backscatter. Since the soil is dry and the radar wavelength is 
long, deep penetration into the soil is expected. Because "Sp2" has the low stand density 
and "St2" has small trees, more soil surface is exposed in these stands. Thus, the back- 
scatter from the sub-surface may be a major contributor. 

The model predicts well for the HH and VV correlation coefficients of stands 
"St2" and "Stll", but overestimates the coefficient of stand "Sp2" (Figure Id). The 
overestimate can be attributed to inadequate modeling of surface backscatter. The sur- 
face model used gives completely correlated backscatter for HH and VV because identi- 
cal surface parameters are used for all the simulated pixels in a stand, and because the 
surface model predicts zero cross-polarized backscatter. For the sparse stand "Sp2", the 
surface backscatter is the major scattering source. Thus, strong and completely corre- 
lated surface backscatter is predicted by the model. 

The model predicts a phase difference of « 210.0° for the dense stands "St2” and 
"Stll". Since the model phase difference and the SAR observed phase differences are 
close to 180° for stands M St2" and "Stl 1", there is some indication that the VV-HH phase 
difference from double-bounce trunk-ground interactions (180°) contributes to the total 
phase difference. 

To investigate the model prediction quantitatively, we test the null hypothesis of 
H a : |i = \x sar > where is the modeled mean and p jar is the mean of the SAR 
data. At the 5% significance level, the H 0 tests of the HH, HV, and VV backscatter, the 
correlation coefficient of HH and VV polarizations, and the VV-HH phase difference are 
accepted. 

5. Conclusions 

For P-band, the model predicts well compared with SAR data for HH, HV, and VV 
backscatter, VV-HH phase difference, and HH and VV correlation coefficient for stands 
with medium-to-high density (228 - 303 trees per hectare). When stand density is low 
(23 trees per hectare), model performance becomes unacceptable: the model overesti- 
mates HH and VV backscatter, underestimates HV backscatter, and overestimates HH 
and VV correlation coefficient. 

6. References: 

[1] J. A. Richards, G. Sun and D. S. Simonctt, fi L-band radar backscatter modeling of 
forest stands", IEEE Trans . on Geosci. and Remote Sensing, vol. GE-25, no. 4, pp. 
487-498, 1987, 

[2] M. A. Karam and A. K, Fung, "Electromagnetic scattering from a layer of finite 
length, randomly oriented, dielectric, circular cylinders over a rough interface 
with application to vegetation". Ini . J. of Remote Sensing, vol. 9, no, 6, pp. 1109- 
1134,1988. 

[3] S. L. Durden, J. J. van Zyl, and H. A, Zcbkcr, "Modeling and observation of the 
radar polarization signature of forests areas", IEEE Trans, on Geosci . and Remote 
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Sensing, vol. 27, no. 3,pp. 290-301, 1989. 

[4] F. T. Ulaby, K. Sarabandi, K. C. McDonald, M. Whitt, and M. C. Dobson, "Michi- 
gan Microwave Canopy Scattering Model (MIMICS) , Int. J. of Remote Sensing, 
vol. 1 1, no. 7, pp. 1223-1253, 1990. 

[5] G. Sun, D. S. Simonett, and A. H. Strahler, "A radar backscattering model for 
discontinuous coniferous forests", IEEE Trans, on Geosci. and Remote Sensing, 
vol. 29, no. 4, pp. 639-650, 1991. 

[6] Y. Wang, G. Sun, and John Day, "Santa Barbara microwave backscattering model 
for woodlands", submitted to Int. J. Remote Sensing for publication, pp. 27, 1991. 

[7] Y. Wang, J. L. Day, F. W. Davis, and J. M. Melack, "Modeling L-band radar back- 
scatter of Alaskan boreal forest", submitted to IEEE Trans, on Geosci. and Remote 
Sensing, for publication, pp. 19. 

[8] Y. Wang, "Radar backscattcr canopy modeling and applications", A Ph. D. Disser- 
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Figure 1 Model results vs SAR data at P-band 



-14 -10 -8 -6 -4 

(a). Model HH backscatier 


-16 -14 -12 -10 -8 

(b). Model W backscatter 



Solid lines indicate the equality or modeled means to those of SAR data. Dotted 
lines show ± 1.0 dB calibration uncertainty for the HH, HV, and VV backscatter (Figure 
la-c) and ±0.1 for the HH and VV correlation coefficient (Figure Id). "+", and "x" 
represent data points of stands "Sp2", "St2", and "Stll", respectively. 
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Introduction 

A study is in progress comparing AIRSAR backscatter from coniferous forest 
plots containing gaps to backscatter from adjacent gap-free plots. We are asking "How 
do gaps in the range of 400 to 1600 m 2 (approximately 4-14 pixels at intermediate 
incidence angles) affect forest backscatter statistics?" and "What incidence angles, 
wavelengths, and polarizations are most sensitive to forest gaps?" In order to visualize 
the slant-range imaging of forest and gaps we make use of a simple conceptual model. 
This strictly qualitative model has led us to hypothesize that forest radar returns at 
short wavelengths (eg., C-band) and large incidence angles (eg., 50°) should be most 
affected by the presence of gaps, whereas returns at long wavelengths and small 
angles should be least affected. Preliminary analysis of 1989 AIRSAR data from 
forest near Mt. Shasta supports the hypothesis. 

Current forest backscatter models such as MIMICS (Ulaby et al. 1990) and 
Santa Barbara Discontinuous Canopy Backscatter Model (Sun et al. 1991) have in 
several cases correctly predicted backscatter from forest stands based on inputs of 
measured or estimated forest parameters (McDonald et al. 1990, Wang and Paris 
1992). These models do not, however, predict within-stand SAR scene texture, or 
"intrinsic scene variability" as Ulaby et al. (1986) has referred to it. For instance, the 
Santa Barbara model, which may be the most spatially coupled of the existing 
models, is not truly spatial. Tree locations within a simulated pixel are distributed 
according to a Poisson process, as they arc in many natural forests, but tree size is 
unrelated to location, which is not the case in nature. Furthermore, since pixels of a 
simulated stand are generated independently in the Santa Barbara model, spatial 
processes larger than one pixel are not modeled. Using a different approach, Oliver 
(1991) modeled scene texture based on an hypothetical forest geometry. His simu- 
lated scenes do not agree well with SAR data, perhaps due to the simple geometric 
model used. 

Insofar as texture is the expression of biological forest processes, such as succes- 
sion and disease, and physical ones, such as fire and wind-throw, it contains useful 
information about the forest, and has value in image interpretation and classification. 
Forest gaps are undoubtedly important contributors to scene variance. By studying the 
localized effects of gaps on forest backscatter, guided by our qualitative model, we 
hope to understand more clearly the manner in which spatial heterogeneities in forests 
produce variations in backscatter, which collectively give rise to scene texture. 

Shasta Forest Gaps: a Conceptual Model 

The forest at the Mt. Shasta study site spans a range of tree sizes and stand den- 
sities, but for the purpose of a conceptual model, we assume crowns are interlocked 
and tree heights are 20-25 m. Forest gaps ranging from 400 to 1600 m 2 are assumed 
circular. Slant range pixel spacing is 6.7 m. 

The imaging geometry is illustrated in Figure 1 for 0 O =45°. One slant range- 
pixel is shown shaded. If canopy penetration is deep, as at P-band, the entire pixel 
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volume can contribute to backscattcr. In this case, effective range resolution becomes 
a function not only of 0 O and slant-range resolution, but also of tree height, canopy 
depth, and extinction coefficient. Effective range resolution for long wavelengths and 
small incidence angles could be much coarser than ground range resolution. Con- 
versely, at large incidence angles the pixel is oriented more vertically, which reduces 
the pixel volume intersecting the canopy. At short wavelengths, canopy extinction 
limits backscatter to the upper canopy, further reducing pixel volume. It follows that 
short wavelengths combined with large incidence angles should lead to the maximum 
effective range resolution of laterally oriented forest features (such as gaps), and to the 
greatest effects on local backscatter. 


Figure 1 



resolution 


Methodology and Preliminary Results 

To test this hypothesis, 238 forest gaps were located on orthophotoquads of the 
region. Gap coordinates in three SAR images 40°, and 50 ) were deter- 

mined by co-registering a digitized map of the forest gaps to the SAR images. We 
are currently comparing backscatter from plots containing gaps in several size classes 
to adjacent areas without gaps. The mean and median of the backscatter distributions 
from 5x7 pixel windows have proven insensitive to the presence of gaps in the 500- 
900 m ^ size range. A more effective test statistic is the difference in backscattcr at 
the lower quartiles of the gap and non-gap backscattcr distributions. Tests on a subset 
of the data (n=12 pairs) indicate that the lower quartile backscattcr at C-band, HH 
polarization, for 0 o =5O° is significantly lower (at the 95% confidence level) for gap- 
containing plots than for gap-free plots (Table 1). Lower quartile differences for other 
bands and angles are not significant. Following analysis of the full data set, we plan 
to test other parameters, such as the difference in the coefficient of variation of gap 
and non-gap windows, so as to build a clearer empirical understanding of the effects 
of forest spatial heterogeneities on backscatter. 
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Table 1. Difference between backscaitcr from gap-containing and gap-free forest plots 


C-band L-band P-band 


HH HV VV 

HH HV VV 

HH HV VV 

25° 

40° 

50° 

0.56 0.65 0.41 

-0.39 -0.09 -0.38 

*-1.42 -0.85 -0.40 


0.27 0.06 -0.22 

-0.45 -0.01 -0.76 

-0.36 0.29 -0.86 


Values shown are the mean differences in backscaitcr (dB) between the lower quartiles 
of gap-containing plots and the lower quartiles of gap-free plots. * indicates 
significance at the 95% confidence level for the 12 pairs of plots tested. 
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1. INTRODUCTION 

The analysis of JPL AIRSAR data over the Landes forest in South-West 
France revealed strong correlation between L- and especially P-band a° and the pine 
forest biomass (Le Toan et al 1992). To explain the physical link of radar backscatter to 
biomass, a polarimetric backscattering model was developed and validated (Hsu et al. 
1992). Then the model was used in a simulation study to predict <f sensitivity to 
undesired canopy and environmental parameters (Beaudoin et al. 1992). This paper 
reports main results concerning the data analysis, modeling and simulation at P-band. 


2. EXPERIMENT AND DATA ANALYSIS 

The test site is a managed forest of maritime pine (Pinus Pinaster), including 
dear-cuts and forest stands from seedlings up to 46 years old, providing a wide range of 
forest parameters (height, dbh, tree density, etc.). The collected parameters concern 
forest inventory, measurements during radar overflight and tree structure sampling. P- 
L-C band quad-pol AIRSAR data were acquired on August 16, 1989. A complete 
polarimetric calibration procedure was appplied to the data. 

Of particular interest for applications are the above-ground biomass by parts 
(needle, branch and trunk biomass), which can be estimated using tree basal area and 
tree density. Fig. 1 shows biomass by parts which is increasing quasi linearly with age. 
High correlation was found between L- and especially P-band <f and most forest 
parameters, especially biomass (tons/ha), for example total above-ground biomass as 
shown in Fig. 2. As <f was found statistically correlated to many forest parameters (age, 
height, dbh, biomass), which are themselves biologically interrelated, theoretical 
modeling is needed to explain the interaction between incident wave and different parts 
of forest canopies. 

3. MODELING 

A radiative transfer model is used for the modeling of the pine forest 
backscatter (Hsu et al. 1992), which includes a branching model for vegetation clusters 
found in pine forests. The forest is modelled as a 4-layer discrete random medium over 
a slightly rough ground (crown, trunk, bush and grass layers). Vegetation constituents 
are modeUed as clusters of dielectric cylinders accounting for needles, twigs, branches 
and trunks. Both incoherent and coherent scattering from these multiscale clusters are 
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computed. The model gives the expression of the lst-order solution to the RT equations 
as the sum of volume and surface-volume scattering terms in each layer. For hilly 
terrains, the model considers the ground slope and aspect angles. Fig.2 shows 
model/measurement comparison of P-band <f at HH,HV and W polarizations at 
0=45° , as a function of the total above-gound biomass, while Fig. 3 presents the main 
scattering contributions at HH and W. The HH return is mainly governed by crown 
and crown-ground scattering for young forest, and trunk-ground scattering for mature 
forest. The W and HV returns are dominated by crown scattering, which comes from 
the primary branches. Thus, W and especially HV o° are directly linked to radius, 
length, density and moisture content of primary branches. Hence o° increases with 
increasing branch fresh biomass. Relations between a° and forest height, trunk or total 
biomass are therefore indirect. The above results provide a significant improvement of 
our knowledge on the scattering mechanisms and consequently our understanding of the 
relationships <f- forest biomass. 

4, CHANGING FOREST AND ENVIRONMENTAL PARAMETERS 

To extend the above observations to various forest and environmental 
conditions, a simulation study is performed. 

1) Ground slope effects: Fig.4 presents <f as a function of local ground slope facing the 
radar. W and HV are not changed significantly, whereas HH is fastly decreasing for 
old forests where the trunk-ground term dominates, and is approaching to the crown 
scattering term. 

2) Varying understorv conditions : <f is computed for a canopy without and with an 
understory vegetation (lm high, LAI = 1.0). As expected, HV and W returns are not 
affected, whereas a ldB decrease of <f for HH is due to the attenuation of the trunk- 
ground term. 

3 Wegetation moisture condition was varied from 55% to 65% to simulate seasonal or 
diurnal effects. HH is not significantly sensitive to this variation, whereas W and HV 
returns are more affected (about ldB). 

3) Crown structure: <f was simulated as a function of 0, the mean elevation angle of the 
primary branches, which is a major species-dependent structural parameter. At HV and 
W, the maximum sensitivity to crown biomass is found in the range 45°-60°. For low or 
high 0 (erectophile or planophile species), there is no significant sensitivity to crown 
biomass. For ? HH returns, no effect of 0 on the a° sensitivity to biomass is observed, as 
the return results mainly from trunk-ground interaction. 

5. DISCUSSION & CONCLUDING REMARKS 

At present, the modeling and simulation studies were performed at P-band, 
which was found optimal for retrieval of forest biomass. HH return was found governed 
by trunk-ground scattering, and thus can be significantly affected by environmental 
conditions such as soil moisture, ground slope and bush layer. On the other hand, W 
and especially HV returns were found tightly linked to crown biomass and thus 
unaffected directly by environmental conditions. Therefore, it is possible to derive total 
above-ground biomass from P-band HV or W SAR data, knowing relations between 
biomass by parts. HH return offers the advantage of being less species dependent. 
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1. Introduction 

The Shuttle Imaging Radar-C/X-SAR Experiment 
includes the study of wetland dynamics in the 
seasonal tropics. In preparation for these wetland 
studies, we analyzed airborne P, L, and C band 
radar (AIRSAR) data of Belize, Guatemala, and 
Mexico acquired by NASA and JPL in March 1990. The 
first phase of our study focuses on AIRSAR data 
from the Gallon Jug test site in northwestern 
Belize, for which ground data were also collected 
during the three days prior to the overflight. One 
of the main objectives of the Gallon Jug study is 
to develop a method for characterizing wetland 
vegetation types and their flooding status with 
roultifrequency polarimetric radar data. 

2. Methods of AIRSAR Analysis 

We converted the four look AIRSAR data 
acquired from JPL in the compressed Stokes matrix 
format into five re-scaled 8-bit rasters for each 
frequency (Table 1) . The image was rectified to 
ground range. The HH and W rasters were corrected 
for scene and sensor effects caused by changes in 
look angle across track. This was accomplished by 
compiling backscatter means for three sets of 7 
sample polygons of homogeneous upland forest 
distributed across track, and then performing a 
regression analysis to determine the angular 
dependency of HH and W. The regression equation 
was then applied to correct the image for variable 
look angles. Next, we filtered each raster with 
two iterations of a 3 x 3 modal filter to reduce 
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the effects of speckle. Since we were primarily 
interested in characterizing vegetation types, we 
transformed the AIRSAR parameters from each band 
into four indices that reflect aspects of the 
vegetation (Table 1) . These indices (except for 
the biomass index) have the added benefit of being 
mostly independent of terrain slope, thus allowing 
for better measures of vegetation characteristics. 
The canopy structure index provides a measure of 
the dominance of vertical structures such as stems 
or trunks. The volume scattering index provides a 
measure of volume scattering, which can be related 
to canopy thickness. The biomass index reflects 
the relative amount of woody compared to leafy 
biomass in forests and therefore has higher values 
for deciduous than for evergreen forests due to 
absorption by green leaves, especially in L and P 
band. The biomass index for herbaceous types more 
closely reflects total above ground biomass. The 
interaction type or phase angle difference index 
measures the importance of single bounce or volume 
scattering relative to double bounce interactions. 
The possible 360 degrees of phase shift have been 
scaled in a monotonic fashion so that interactions 
containing appreciable amounts of double bounce 
have lower values. It is important to note that 
these indices must be interpreted with caution, 
especially in comparisons between ecosystem types 
with very different structure. 

The final step was a clustering analysis. 
Clustering was first run on a pixel by pixel basis 
for the four indices from P, L, and C band using 
the Isoclass algorithm with 45 classes. A second 
clustering of the 45 classes, based on class means, 
was performed using the centroid method in the SAS 
statistical package. 

3 . Results 


The cluster analysis identified 12 distinct 
land cover types in the Gallon Jug test site that 
were verified with field data: 1) pasture and 
cultivated lands, 2) water, 3-4) two types of 
marsh, 5) recently cut forest, 6) regrowth, 7) 
flooded swamp forest, 8) non-flooded swamp forest, 
9) flooded swamp thicket, and 10-12) three types of 
upland forest. Examination of the indices for 
each of the 12 types helps clarify differences. 
The two marsh types represent low (C band BM=29) 
and high (C band BM=92) biomass marshes. The swamp 
forests have a high percentage of woody compared to 
leafy biomass (P band BM=63-72) , especially when 
swamp forests are compared with upland forests (see 
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below) . Swamp forest flooding is indicated by the 
lower PD (P band PD=43) as compared with the non- 
flooded swamp forest (P band PD=62) . The swamp 
thicket has less woody biomass (P band BM=33) than 
the swamp forests, and again flooding is indicated 
by the low PD (P band PD=45) . The three forest 
types represent a gradient from semi-deciduous to 
semi— evergreen as indicated by decreasing woody 
biomass (P band BM =46, 40, 27) . This decrease in 
BM is mirrored by an increase in P and L band 
volume scattering (e.g. P band VS = 49, 50, 52), 
which may reflect an increase in canopy closure or 
thickness. The BM index for the agricultural 
ecosystems shows an increase from pasture/ 
cultivation (e.g. P band BM=6) to regrowth (e.g. P 
band BM— 14) , although some of the dense forest 
regrowth was classified as forest. The recently 
cut forest has felled trees with no leaves, which 
give a very high woody biomass value (P band 
BM=85) , a low volume scattering value (P band 
VS=40) , and a low canopy structure index (P band 
CAN=40) . 


Table 1. AIRSAR parameters and indices used in 
analysis (compiled for P, L, and C band) . 

Initial parameters derived from Stokes matrix* 

Backscatter amplitude parameters 

HH W CS = HV+VH LK = W + HH 

2 2 

Backscatter phase angle parameter 
HH,W phase angle difference = PD 

AIRSAR indices 


Canopy structure = CAN 
Volume scattering = VS 
Biomass = BM 


W 

W + HH 
CS 

CS + LK 
LK 


Interaction type = PD = pd converted to 
(HH,W phase angle diff.) roonotonic scale 


ik * t 

Derivation of these parameters is described in 
detail in J . F. Paris, Processing Airborne Imaging 
Radar (AIR) Data on a Microcomputer Workstation. 
Technical Report TR-89-04 , Department of Geography, 
California State University, Fresno, 1989. 
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Introduction 

Mapping of floodplain inundation patterns is a key element in developing hydro- 
logical and biogeochcmical models for large tropical river basins such as the Amazon. 
Knowledge of the time sequence of inundation is necessary to determine both water 
routing and biogenic gas fluxes (Richey et al. 1990). SAR is uniquely su.lcd for this 
application because of its ability to penetrate cloud cover and, in many cases, to detect 
flooding beneath a forest or herbaceous canopy (Hess et al. 1990). We are currently 
developing a procedure for discriminating flooded forest, flooded herbaceous vegetation, 
and open water from other cover types for a coastal wetland site on the lower Altamaha 
floodplain, Georgia, emphasizing robust classifiers that are not site-specific. 

Study site and methods 

Multiple dalatakes over a range of incidence angles were obtained for the 
Altamaha site by the JPL polarimclric SAR in March 1990 and again in May 1991. For 
both dates, the Altamaha was at high flood stage, and the entire floodplain other than 
sand ridges was inundated. Water levels were documented by extensive ground observa- 
tions during the overflight. A rich variety of wetland types occur m the study site, 
including cypress-tupclo swamp (Taxodium distichum. Nyssa aquaiica). bottomland 
hardwood forest (Nyssa spp., Taxodium distichum, Fraxinus spp., Acer rubrum, Liqui- 
dambar styraciflua), and marshes dominated by Spartina allerniflora, Juncus roemeri- 
anus, and Zizaniopsis miliacea. Upland sites are predominantly loblolly pine (/W 
laeda) plantations. Results given in this summary are fora March 1990 scene calibrated 
with JPL’s SUPERCAL software using comer reflectors. The radar variables used in the 
analysis are HH, VV, and HV radar cross-sections per unit area (a 0 ) and HH-VV phase 
difference for C, L, and P-bands, and ratios directly derived from these parameters. A 
5X5 median filter was applied to the radar cross-sections to reduce within-class variabil- 
ity. Initial analysis was performed on 100 pixels in each of five cover types: open water, 
clearing, marsh, pine, and flooded forest. For each category, five 20-pixel windows were 
selected, representing a range of incidence angles and vegetation variability within cac 
category. Because nearly all hardwood forests in the scene were flooded at the time ot 
the overflight, pine forests are used to represent unflooded forests. 


Discrimination of cover types 

Discussion in this summary will be limited to two of the radar parameters found to 
be the most useful: PHH to CHH ratio and L-band phase difference. In Figure \,a P „„ is 
plotted vs. aam for unfiltcred and median-filtered data; each plotted character 
represents a pixel. Separation between the classes is good, though not absolute, for the 
filtered data- overlap in the unfiltered data makes discrimination impossible. For the 
filtered data, mean o° values for flooded forest and pine respectively are -8.2 dB and 
-12 8 dB at CHH and -1.1 dB and -4.6 dB at PHH. Although the means are distinct at 
both bands, there is enough spread in the data to cause confusion between the pine and 


21 


flooded forest categories. 

Figure 2 demonstrates the behavior of phase difference at L band and C band for 
flooded forest. These phase difference histograms are presented on a circular scale 
rather than the usual linear one. There is a distinct shift toward phase differences close to 
180° for flooded forests at L band, indicative of double-bounce scattering; this is not the 
case for pines. However, the phase difference distribution is too broad to perform well 
as a classifier, and phase difference, because of its circular nature, does not lend itself to 
techniques such as median filtering. We have found that by classifying the phase 
difference into one of three types and then applying a modal filter, confusion between 
pines and flooded forests is largely eliminated. Phase differences from 0°-70° or 270°- 
360° are classed as type 1 , those centered on 180° (1 10°-250°) are considered type 3, and 
intermediate phase differences (70°-110° or 250°-290°) are type 2. A 3X3 pixel modal 
filter is then applied to the data, creating regions of uniform phase difference type. 
Although this technique is conceptually very simple, results to date show excellent 
separation between flooded forest and other categories. Table 1 shows the number of pix- 
els classified in each type for 1000 pixels each from the marsh, pine, and flooded forest 
categories. Ninety-nine percent of the flooded forest pixels were classed as type 3, Only 
1 non-flooded-forest pixel was type 3, and six flooded forest pixels were type 2. 


Table 1. LHH-LVV phase difference types for 1000-pixcl 

samples of marsh, pine, and flooded forest 


Type 1 

Type 2 

Type 3 

Marsh 

520 

479 

1 

Pine 

970 

30 

0 

Flooded forest 

0 

6 

994 


The scatter plot for filtered Cp/ m vs. Genu suggests that, following widely used 
parametric classification procedures, combinations of parameters could be found to pro- 
vide a highly accurate separation between the cover types. The classification procedure 
would be tailored to the radar signatures (means and covariance matrices) for that partic- 
ular scene. The limitation of this approach is that good results usually would not be 
obtained when applying the signatures to other scenes. Even for the same scene, signa- 
tures will vary temporally due to seasonal changes in scene properties such as leaf area 
index. Because we plan to apply this procedure to a broad spatial and temporal domain, 
scene-by-scene optimization of parametric classifiers is impractical. A knowledge-based 
approach using convergent evidence (Wharton 1989) seems best suited to this task. 
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Fig. la. Unfiltered data. Clearing (C), Marsh (M), Pine (P), Flooded Forest (F) 



Fig. 1b. Filtered data, Clearing (C), Marsh (M), Pine (P), Flooded Forest (F) 



Fig. 2a. HH-W phase difference Fig. 2b. HH-W phase difference 
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L INTRODUCTION 

We present a study on the performance of a supervised fully polarimctric maximum 
likelihood classifier for S AR data when applied to a specific classification context: for- 
est classification based on age classes and in the presence of a sloping terrain. 

For the experimental part of the study we use the polarimctric AIRS AR data at P, L, 
and C band, acquired over the German Black Forest near Freiburg in the frame of the 
1989 MAESTRO-1 campaign and the 1991 MAC Europe campaign. MAESTRO-1 
was an ES A/JRC sponsored campaign, and MAC Europe (Multi-sensor Aircraft Cam- 
paign) was organized by NASA; in both cases the multi-frequency polarimetric JPL- 
AIRSAR radar was flown over a number of European test sites. 

The study is structured as follows. At first we investigate the general characteristics of 
the classifier and the dependencies from some parameters, like frequency bands, fea- 
ture vector, calibration, using test areas lying on a fiat terrain. Once we have determined 
the optimal conditions for the classifier performance, we then move on to the study of 
the slope effect. The bulk of this work is performed using the Macstrol data set. 

Next we consider the classifier performance with the MAC Europe data. The study here 
is divided into two stages: first we repeat some of the tests done on the Maestro data, 
to highlight the improvements due to the new processing scheme that delivers 16 look 
data. Second we experiment with multi images classification with two goals: to assess 
the possibility of using a training set measured from one image to classify areas in dif- 
ferent images; to classify areas on critical slopes using different viewing angles. 

In this summary paper we will list the main points of the study and highlight only some 
of the results obtained so far. 

2. CLASSIFIER PRINCIPLES 

The principles of the fully polarimetric maximum likelihood classifier are described 
in [Kong et al . 1 987] . S uffice here only to say that the discrim inant function, in the form 
of a minimum distance, is given by: 

djx) = x + • C* 1 • x + In I CJ, (1) 

where x = j HH HV wj is the complex feature vector for the unknown pixel and C m 
is the spatially averaged covariance matrix of the training class. 

3. GENERAL CLASSIFIER CHARACTERIZATION 
3.1 Methodology 

For the study of the classifier performance, we define training and test sets using the 
ground data collected during the Maestrol campaign [Churchill 1989]; the forest is 
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segmented into stands characterized by the mean tree age. The tree ages have been 
grouped into four main classes: young trees (10-30 years), mature trees (30-70 years), 
adult trees (70-110), old trees (>110 years).The training and test areas are taken on a 
flat terrain, to filter out from this analysis the influence of slope. 

At this level of the study we then investigate the following points (see also [De Grandi 
etal. 1992]): 

• characterization of the training sets through their covarian- 

ce matrices 


• classification performance as a function of the frequency 
bands 

• dependency from the off diagonal terms of the covariance 
matrix 


• dependency from polarimetric calibration 

• relative weight of the feature vector elements in the discri- 
minant function 

• performance in the case of a one channel and a two chan- 
nels feature vector 


3.2 Results and discussion 

As an example we report in table I the results of the classification of 3 forest age classes, 
an urban area and an arable field at P, L and C bands. 

Table I - Confusion Matrix for 3 Forest Age Classes, an Arable Field and an 
Urban Area; Maestrol Freiburg Data Set at P-, L- and C Bands 


r 

— — : 1 


Pband 

L band 




L band 




urban 1 

arable 1 

agd 1 

ag e 2 1 

age4 




age2 


urban 

arable 

agel 

age2 

agc4 

BSS 











0.0 


9.5 

9.2 

2.6 
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■u 

9.3 
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■9 

0.0 

14.6 

34.2 

w 

29.3 






wm 


■t|| 




0.0 


21.1 

K m 

30.1 

age4 







1123 

KaU 

28.3 

i 

0.0 

■S 

19.7 

| 26.1 

46.8 


Given the classification parameters of our test case, P band appears to be the most suit- 
able for forest discrimination. 


4. SLOPE EFFECT 


4.1 Methodology 

The Maestrol Freiburg data set has been georeferenced [De Groof ct al. 1991] and it 
is therefore, possible to attach to each pixel in the image space the corresponding coor- 
dinates in the map space, and additional layers of information, such as the local inci- 
dence angle, and the slope. This processing step has then allowed us to divide each 
polygon representing an age class into sets of pixels having a given slope angle. In par- 
ticular we defined 6 slope classes with angular spacing of 5 degrees. This spacing was 
chosen from considerations on the variance of the powers of the 3 polanmetnc chan- 
nels. 


4.2 Discussion 

The test set for 3 age classes exhibits high between groups variance when the 6 slope 
classes are considered and the confusion in the classification is therefore relevant. A 


25 









































possible correction to the slope influence could be to drive the classifier with the knowl- 
edge of the slope angle for each pixel; in other words, the classifier could be switched 
among a number of different training sets, according to the slope angle of the pixel to 
be classified. Unfortunately this method works only for certain slope classes; indeed 
for positive slope angles greater than 10 degrees, the class separation is also affected, 
and therefore the classification accuracy is degraded even if a proper training set is tak- 
en into account. This can be seen for instance calculating the distance between 2 age 
classes and the 6 slope classes. 

5. MAC EUROPE DATA 

5.1 Methodology 

The MAC Europe data have been processed with a new correlator and one of the avail- 
able products is a 1 6 look data set Moreover for the Freiburg test site we have at dispos- 
al a number of flights paths, both in the same direction and at orthogonal directions. 
This paves the way to an interesting experiment in multi image classification. First and 
essential step for this purpose is an accurate radiometric and polarimetric calibration 
of the data sets, and their co registration. By the time of writing of this summary this 
processing step is under way. 

The plans are to explore the possibility of classifying multiple scenes using a training 
set measured on one scene, and to exploit the different viewing angles to resolve the 
confusion in forest age classification introduced by the terrain slope. 

5.2 Preliminary results 

Classification based on the same training set used in the test case with the Maestro data 
reported in section 3.2 has been performed on the Mac Europe data set with scene iden- 
tifier cm3207. The results are reported in table II. The classification accuracy is signifi- 
cantly improved at P band with respect to the Maestrol case. 

Table II - Confusion Matrix for 3 Forest Age Classes and an Arable Field; 
MAC Europe Freiburg Data Set at P-, L- and C Bands 



Pband 

Lband 

C band 


arable 

agel 

age2 

age4 

arable 

agel 

age2 

age4 

arable 

agel 

age2 

age4 

arable 

96.69 

1.24 

0.00 

0.00 

nmol 

2.89 

Q | 

0.41 

92.15 

3.72 

0.41 

2.07 

agel 

0.00 

92.86 

3.57 

3.57 

EESI 

64.29 

Egj 

3.57 

21.43 

57.14 

17.86 

3.57 

age2 

0.00 

9.90 

86.46 

3.65 

0.00 

38.02 

54.69 

7.29 

25.52 

7.81 

42.19 

24.48 

age4 

0.00 

1.06 

8.46 

90.48 

0.00 

56.34 

28.85 

14.80 

34.74 

14.50 

13.29 

37.46 
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1. INTRODUCTION 

During the last several years, significant efforts in microwave remote sensing 
have been devoted to relating forest parameters to radar backscattering coefficients 
(e.g., Le Toan et al., 1991; Kasischke, et al., 1991). These and other studies showed 
that in most cases, the longer wavelength (i.e. P band) and cross-polarization (HV) 
backscattering had higher sensitivity and better correlation to forest biomass. 

This research examines this relationship in a northern forest area through both 
backscatter modeling and SAR data analysis. The field measurements (Ranson and 
Smith 1990) were used to estimate stand biomass from forest weight tables. The 
backscatter model described by Sun et al. (1991) was modified to simulate the back- 
scattering coefficients with respect to stand biomass. The average number of trees per 
square meter or radar resolution cell, and the average tree height or diameter breast 
height (dbh) in the forest stand are the driving parameters of the model. The rest of 
parameters such as the dielectric constants of tree components and soils, roughness of 
the soil surface, orientation and size distributions of leaves and branches, remain 
unchanged in the simulations. 

2. MODEL DESCRIPTION 
2.1. Tree crown scattering 

Tree crowns were modeled as ellipsoids consisting of a mixture of leaves and 
branches (scattering elements) with various sizes and orientations. The scattering can 
be characterized by the complex scattering matrix S. The Stokes matrix L derived 
from the complex scattering matrix S was averaged over all orientations and sizes of 
the scattering elements. 

The extinction coefficient matrix k is expressed in terms of the forward scatter- 
ing amplitude, and was averaged over orientation and size distributions of scatterers in 
the same way as the Stokes matrix. Having multiplied by the number of scatterers per 
cubic meter in the tree crown, the Stokes and Extinction matrices then can be used in 
the radiative transfer equation to calculate the backscattering and attenuation from a 
tree crown layer. 

Table 1 and 2 show some calculated results using input data from a hemlock 
forest stand (Chauhan et al., 1991). The first element of Stokes matrix 
(L(l,l) = S vv S V v*) was used lo compare the backscattering strength between leaves 
and branches, and between branches with various sizes. The first element of extinc- 
tion matrix (k( 1,1) = 2Re(z XS vv f ) was used to compare attenuation from various com- 
ponents within the tree crown. 

Table l shows that a) attenuation by leaves is higher than branches in all bands; 
b) backscattering by branches is higher than leaves except at C band; c) backscattering 
by leaves at P (and L) band is negligible; and d) to account for scattering and 



attenuation both leaves and branches must be considered in the model. 

Table 2 shows that backscattering of branches at C band is mostly contributed 
by secondary branches, i.e. AA, BB and CC; at L band by A, B, and C branches, and 
at P band by the large primary branches, i.e. B, C and D branches. 

Table 1. Backscattering and attenuation of leaves and branches. 196000 needles and 
105.3 branches per m 3 were assumed. 


Type 

C Band 

L.I K„ 

L Band 

Lb *u 1 

P Band 

I'll K n 

Leaves 

2.39e-3 

0.6384 

1.02e-5 

0.13348 

2.31e-7 

0.07630 

Branches 

1.47e-3 

0.1158 

1.0e-3 

0.06862 

5.57e-4 

0.003876 

Table 2. 

— 1 

Branch physical parameters and their relative volume backscattering (vv) 

strength. 

Percentages reflect branch size probabilities. 



Type 


Diameter 

length 

Proba- 


Relative <7 (%) 




(cm) 

(cm) 

bility 

C Band 

L Band 

P Band 

AA 


0.2 

15 

0.86684 

13,66 

0.26 

0.0 

BB 


0.4 

30 

0.10312 

43.49 

1.26 

0.09 

CC 


0.6 

40 

0.01605 

14.88 

1.47 

0.11 

DD 


0.8 

75 

0.00619 

8.01 

3.88 

0.38 

A 


1.1 

90 

0.00481 

5.15 

14.05 

1.3 

B 


2.4 

165 

0.00154 

4.65 

35.64 

18.90 

C 


3.0 

260 

0.00119 

7.44 

37.42 

54.81 

D 



3.5 

350 

0.00025 

2.72 

6.03 

24.57 


2.2. Model simulation 

The following major backscattering components were considered by the model: 
1) direct backscattering from tree crowns; 2) direct backscattering from ground surface 
attenuated by tree canopy; 3) crown-ground double-bounce scattering; 4) trunk-ground 
double-bounce scattering; and 5) direct backscattering from trunks. 

In the model the ground surface was characterized as very rough - with the stan- 
dard deviation of the surface height (rms height) of 24.3 cm, and correlation length of 
1.0 meter. A geometrical optics model was used to calculate both the backscattering 
and the scattering in the specular direction. Because of the high attenuation from the 
tree canopy and the roughness of the surface, the components 3 and 4 were very low, 
even in P band. The dominant component is the direct backscattering from the tree 
crown, with some contribution from component 2 at HH and VV polarizations, espe- 
cially in low biomass and long wavelength cases. The HV backscattering was solely 
from tree crown scattering. A random leaning angle (normal distribution with mean 
of zero and standard deviation of 5 degrees) was assigned to each tree, but the contri- 
bution to HV backscattering from the direct trunk backscattering was not significant. 

Figure 1 shows some simulation results. As mentioned earlier, the number of 
trees per pixel and average height or dbh are inputs to backscatter model. Thus, both 
the backscattering coefficients and the total biomass of this pixel can be calculated. 
The lines in these figures show the simulations results. Biomass of several stands were 
calculated from available ground truth and are plotted in the figures along with back- 
scattering coefficients extracted from JPL AIRSAR data of June 9, 1991. 

3. CONCLUSIONS 

Generally, the simulations give satisfactory results, except that the backscatlcr- 
ing from L band is about 2 dB higher than the JPL SAR data. The possible coupling 
between branch sizes and L band wavelength may be the cause of this problem. The 


28 





eight distinct categories of branches are far from representing reality. 

Also shown in the figure are the ratios of P/C and L/C. In another study (Ranson 
and Sun, 1992), it was found that the ratio of HV backscattering from a longer 
wavelength (P or L) to that from a shorter wavelength (C) was a better combination 
for mapping forest biomass. This ratio enhanced the correlation of the image signa- 
ture to the standing biomass and compensated for a major part of the variations in 
backscattering attributed to radar incidence angle. 



Figure 1. The relationships between backscattering and forest biomass: lines indicate 
modeling results, and symbols indicate JPL AIRSAR data. 


4. REFERENCES 

N. S. Chauhan, R. H. Lang, and K. J. Ranson, Radar modeling of a Boreal forest, 
IEEE Transactions on Geoscience and Remote Sensing, Vol. 29, No. 4, 1991, 
pp. 627 - 638. 

E. S. Kasischke, N. L. Cristensen, Jr., and M. C. Dobson, The relationship between 
aboveground biomass and radar backscatter as observed on Airborne SAR 
imagery. Proceedings of the Third Airborne Synthetic Aperture Radar (AIRSAR) 
Workshop, May 23 and 24, 1991, Jet Propulsion Laboratory, Pasadena, Califor- 
nia, pp. 11-21. 

T. Le Toan, A. Beaudoin, J. Riom, and D. Guyon, Relating forest parameters to SAR 
data, IGARSS'91, 1991, pp. 689-692. 

K. J. Ranson, and J. A. Smith, Airborne SAR experiment for forest ecosystems 
research - Maine 1989 Experiment, IGARSS' 90, The University of Maryland, 
College Park, Maryland, 1990, pp. 861-864. 

K. J. Ranson, and G. Sun, Mapping biomass for a northern forest ecosystem using 
multi-frequency SAR data, Proceedings of IGARSS'92, May 1992, Houston, 
Texas, USA. 

G. Sun, D. S. Simonett, and A. H. Strahler, A radar backscatter model for discontinu- 
ous coniferous forests, IEEE Transactions on Geoscience and Remote Sensing, 
Vol. 29, No. 4, 1991, pp. 639-650. 


29 






N94-16742 tj 


SAR OBSERVATIONS IN THE GULF OF MEXICO 

David Sheres 

Center for Marine Science 
University of S. Mississippi 
Stennis Space Center, Mississippi 39529 

1 . INTRODUCTION 

The Gulf of Mexico (GOM) exhibits a wealth of 
energetic ocean features; they include the Loop 
Current with velocities of about 2 m/s and strong 
shear fronts, mesoscale eddies, double vortices, 
internal waves and the outflow of the "Mighty 
Mississippi" river. These energetic features can 
have a strong impact on the economies of the states 
surrounding the Gulf. Large fisheries, oil and gas 
production as well as pollution transport are 
relevant issues. These circulation features in the 
Gulf are invisible to conventional IR and visible 
satellite imagery during the Summer months due to 
cloud cover and uniform surface temperatures. SAR 
imagery of the Gulf does penetrate the cloud cover 
and shows a rich assembly of features there year- 
round. Below are preliminary results from GOM SAR 
imagery taken by SEASAT in 1978 and by the AIRSAR 
program in 1991. 

2. SEASAT IMAGERY 

A search at JPL archive for SAR data 
(optically processed) along the edges of the Loop 
Current yielded images with a rich array of 
features. A number of these images were selected 
and digitally processed. Here we will look at some 
of them with an emphasis on circular doughnut like 
features and their comparison with Shuttle 
photography of what appears to be the same 
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features. These features have a relatively dark 
center surrounded by a bright "doughnut . Their 
spatial scales vary but are generally of the order 
of 10 km. They have been observed before (Fu ana 
Holt, SEASAT Views Oceans and Sea Ice With SAR, JPL 
publication 81-120, 1982) and described as local 

atmospheric downdrafts with the dark center as a 
tropical rain cell. Indeed recent in-sft u 
measurements by Trizna at NRL show the inhibiting 
effect of rain on radar returns confirming the 
above interpretations of these features. 

Large areas in these images are relatively 
bright or dark suggesting large areas with either 
high or low radar crossections. The boundary 
between these areas is a front, most likely the 
edge of the Loop Current. These SEASAT images show 
the power of SAR to monitor energetic features in 
the GOM during the Summer months when other remote 
sensors are practically blind. Thus the modulation 
of waves is a useful means for observing energetic 
ocean features under all-weather conditions. 

AIRSAR IMAGERY OF THE GULF OF MEXICO 

During May 19,1991 we obtained AIRSAR imagery 
from the JPL/NASA DC-8 in conjunction with in-sit a 
data taken from the RV Tommy Munro and tracked 
drifter data. A small number of SAR images were 
processed by the synoptic processor; they do, 
however, show some interesting characteristics. 
During the period of the experiment high winds and 
waves were prevalent. The DC-8 flew a box pattern, 
off the shelf over the Northern GOM, that included 
two equivalent passes with a 45 minute 
difference between them. These passes were to the 
West and in one of them a linear bright front was 
clearly seen while in the other the front was 
missing. The presentation will include wind data 
from the RV Munro that was in the same area at the 
time plus a discussion of different mechanisms for 
the generation of a linear front. This work 
still in progress and I will present as much data 
as will be available at the time of the 
presentation . 
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1. INTRODUCTION 

Extensive AIRSAR measurements were made on 20 July 1990 during the NRL 
Gulf Stream (GS) experiment which addressed a number of scientific questions 
relating to SAR imaging of the ocean surface in the presence of variable currents and 
the background thermohaline circulation (Valenzuela et ah 1991). 

In this paper we concentrate on the nature of the electromagnetic (e.m.) 
backscatter from the North edge of the GS using polarimetric signatures and amplitude 
imagery from the AIRSAR data. 

2. STUDY SITE 

On 20 July 1990 the AIRSAR operated from 9 h 49 min to 15 h 45 min EDT. It 
provided extensive coverage from Norfolk, VA to 37°N, 72°W on the North edge of 
the GS and South across the GS towards the center of a cold ring (34°18'N, 

71°20'W), and return back to Norfolk. For details on the tracks and specific patterns 
performed see Kobrick (1990). 

During the AIRSAR data collection the R/V Cape Henlopen was stationed near the 
North edge of the GS and remained within 30 km of 36.75°N, 72°W the location of 
the main thermal front. The winds ranged from 6 to 10 m/s from the south-west 
(230°T) and waves were almost 1 m high. 

3. POLARIMETRIC SIGNATURES 

The AIRSAR imagery contained a strong linear feature about 20 km south of the 
North edge (thermal front) of the GS and weaker signatures were evident on the South 
edge (Valenzuela et al. 1991). These signatures in the amplitude imagery are strongly 
frequency and polarization dependent. Therefore, a full polarimetric study was done 
to obtain further insight on the nature of the e.m. backscatter for these features. 

Van Zyl et al. (1987) discussed the theory and application of co-polarized 
polarimetric signatures to identify the scattering mechanisms from earth surfaces using 
the shape of the three-dimensional display (Intensity as a function of Orientation and 
Ellipticity angles) as derived from the transmit and receive polarization vectors 
operating on the Stokes matrix of the AIRSAR imagery. 
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In preliminary analysis of polarimetric co-polarized signatures for the North edge 
of the GS (using one frame from pass GS NS 135-1) (Figure 1) and comparisons with 
numerical simulations from a second order tilted-Bragg model (Schuler et al. 1992) 
we have been able to identify that P- and L-band backscatter for the background ocean 
were consistent with Bragg scattering, while for the linear feature and C band other 
contributions are present (Figure 2). 

4. THERMAL FRONTS SIGNATURES 

In addition to the strong linear feature discussed in the previous section, weaker 
signatures were present in the SAR imagery right on the North edge of the GS (the 
main thermal front) and another signature on a thermal front about 7 km north. 

We have extensive sea truth on these fronts for a quantitative correlation with the 
SAR imagery (one frame from pass GS SS 360-1) (Figure 3), but here we investigate 
the frequency and polarization dependence of these signatures for the amplitude 
imagery. From study of the SAR imagery for these features at P, L and C bands and 
various polarizations, it is clear that the signatures are largest for HH and HV 
polarization, and nearly absent for VV polarization. 

Furthermore, we note that the mean intensities for P and L bands do not change 
with azimuth in the image, while C band intensity increases by 2-3 dB towards the 
warm GS waters (Figure 4). The main thermal front is about pixel 800 in the 
azimuthal cut and the other front 7 km to the North is about pixel 200. This increase 
in mean intensity towards the GS is consistent with Ku band scatterometer data from 
the Henlopen and the NRL/P-3 aircraft. 

5. SUMMARY AND CONCLUSIONS 

In this work we have used polarimetric analysis and amplitude imagery to 
investigate the nature of the e.m. backscatter from AIRSAR data near the North edge 
of th GS during 20 July 1990. The e.m. backscatter for the background ocean is 
tilted-Bragg for P and L bands, while C band backscatter contains other mechanisms 
as well, probably specular scattering. 
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Figure 1. JPL ATRSAR total-power image for a 
strong linear ocean feature about 20 km south of the 
North edge of the GS seemingly related to a region of 
large current gradient. 
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Figure 2. A (linear-pol) slice (40° incidence) of the 
3-D polarimetric display of intensity for both the 
background ocean (solid line) and the strong linear 
feature (broken line) of Figure 1. Note the increase 
in horizontal polarization at 0° and 180° angles. 



Figure 3, JPL AIRSAR total-power image for the 
North edge of the GS. The two signatures are 
associated with the temperature fronts, the one on the 
r.h.s. is the main front, the North edge of the GS. 
North is towards the left. 
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Figure 4. Azimuthal intensity cuts (averages of 50 
pixels) of AIRSAR image of Figure 3 for C band and 
various polarizations at about 55° incidence. The 
North edge of the GS is about pixel 800, the other 
temperature front to the north is about pixel 200. 
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1. Introduction 

Under suitable conditions the bottom topography of shallow seas is visible in remote 
sensing radar imagery. Two experiments have been performed to establish which 
remote sensing technique or combination yields optimal imaging of bottom topography 
and which hydro-meteorological conditions are favourable. A further goal is to gain 
experience with these techniques. 

Two experiments were performed over an area in the North Sea near the measuring 
platform Meetpost Noordwijk (MPN). The bottom topography in the test area is 
dominated by sand waves. The crests of the sand waves are perpendicular to the coast 
line and the dominating (tidal-)current direction. A 4x4 km 2 wide section of the test area 
was studied in more detail. 

The first experiment was undertaken on August 16, 1989. During the experiment the 
following remote sensing instruments were used; Landsat-Thematic Mapper, and 
NASA/JPL Airborne Imaging Radar (AIR). The hydro-meteorological conditions: 
current, wind, wave, and air and water temperature were monitored by MPN, a ship of 
Rijkswaterstaat (the OCTANS), and a pitch-and-roll WAVEC-buoy. 

The second experiment took place on July 12, 1991. During this experiment data were 
collected with the NASA/JPL polarimetric SAR, and a five-band helicopter-borne 
scatterometer. Again the hydro-meteorological conditions were monitored at MPN and 
the OCTANS. Furthermore, interferometric radar data have been collected. 

2. Imaging mechanism and implementation 

The imaging mechanism of both the optical and the microwave system is generally 
agreed to consist of the three following steps: 

1. Interaction between bathymetry and (tidal-)current causes modulations in the 
surface current velocity. 

2. Modulations in the surface current give rise to variations in the sea surface 
roughness or wave spectrum. 

3. Variations in the sea surface roughness show up as local modulations in the 
reflected sun light or radar backscatter. 

Based on the above three stage imaging mechanism, an one-dimensional model suite has 
been implemented: 

1 . The interaction between bottom topography and current has been modeled by a two 
stage mechanism. First the depth averaged current is calculated, then the surface 
current is determined from the current profile (Davies 1987 , Zitrnan 1992). 

2. The modulation of the spectrum of the short surface waves through wind, 
dissipation and wave-wave interaction is described by the action balance equation 
(Hasselmann 1960, Willebrand 1975). 

3. Finally, the interaction between sea surface and radar is modelled by a two-scale 
radar backscatter model (Donelan and Pierson 1987, Calkoen et al. 1990). 

3. Data assimilation 

Using the above suite of models, the bottom topography and additional information 
about sea surface temperature, wind speed, and average current, the radar backscatter 
can be predicted. In this paper, interest is focused on the inverse: estimation of the 
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depth from radar backscatter imagery. In principle this can be achieved by explicit 
inversion of the equations. However, this is not feasible, unless some simplifications are 
made, because the equations are non-linear, and noise added to the radar backscatter 
may result in signals outside the range of values attainable by the model. Furthermore, 
additional information obtained from surveys and charts cannot be used. The data 
assimilation approach doesn’t have these disadvantages. It provides the most likely 
bottom topography given all available information. The data assimilation approach 
selects the bottom, which minimizes the difference between measurements and model 
predictions. 

3. Results 

The measurements have been compared with results obtained from the numerical 
models. 

1. The measured depth averaged current and the current profile are in good agreement 
with model predictions (see figure 1). 

2. Sea bottom topography was observed in both optical and microwave imagery. 
However, the practical use of optical data is severely limited by the required sun 
elevation angle and absence of clouds. Sea bottom topography was not observed in 
the C-hand SAR imagery. The usefulness of the recorded L-band imagery is limited 
by dark broad bands caused by interference problems in the radar system. Both HH 
and VV polarized P-band imagery show the bottom topography distinctly. The VV 
polarized P-band imagery corrected for incidence angle dependence is shown in 
figure 2. 

3. Inversion of the model suite by means of the data-assimilation approach is possible. 
The estimated depth based on VV polarized P-band imagery is shown in figure 3. 

4. Conclusions and discussion 

The best depth estimates were obtained from VV polarized P-band imagery. The 
position of the crests of the sand waves can be assessed within 25 metres (two pixels), 
and the height of the crests within 1 m. The errors can partially be explained by the 
time difference between the depth information (1984) and the radar data (1989). 

At present the data obtained during the second experiment are still being calibrated 
and/or pre-processed. The interferometric radar data can be used to determine the 
movements of the water surface. Therefore, it may provide additional information 
useful in the data-assimilation scheme. 
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It is well known that under favorable meteorological and hydrodynamical con- 
ditions the bottom topography of shallow seas can be mapped with airborne 
or spaceborne imaging radar. This phenomenon has been observed for the first 
time in 1969 by de Loor and co-workers in Q-band Side Looking Airborne Radar 
(SLAR) imagery of sandwaves in the North Sea (de Loor and Brunsvelt van Hul- 
ten ’78, Alpers and Hennings ’84, Shuchman et al ’85). 

It is now generally accepted that the imaging mechanism consists of three steps: 

1. Interaction between (tidal) current and bottom topography causes spatial 
modulations in the surface current velocity. 

2. Modulations in the surface current velocity give rise to variations in the 
spectrum of wind-generated waves, as described by the action balance 
equation. 

3. Variations in the wave spectrum show up as intensity modulations in radar 
imagery. 


In order to predict radar backscatter modulations caused by sandwaves, an 
imaging model, covering the three steps, has been developped by the Dutch 
Sea Bottom Topography Group (Vogelzang et al ’89a, Vogelzang et al ’89b, 
Vogelzang et al ’91a). This model and some model results will be shown. 

On August 16, 1989 an experiment has been performed with the polarimetric P-, 
L-, and C-band SAR of NASA/JPL. One scene has been recorded in SAR mode 
(Vogelzang et al ’91b). On July 12 1991 another three scenes were recorded, of 
which one was in the ATI-mode (Along- Track Interferometer). 

These experiments took place in the test area of the Sea Bottom Topography 
Group, 30 km off the Dutch coast, where the bottom topography is dominated 
by sand waves. In-situ data were gathered by a ship in the test area and on 
“Measuring Platform Noordwijk”, 20 km from the centre of the test area. 

The radar images made during the experiment were compared with digitized 
maps of the bottom. Furthermore, the profiles of radar backscatter modula- 
tion were compared with the results of the model. During the workshop some 
preliminary results of the ATI measurements will be shown. 
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1. INTRODUCTION 

The valleys of the Mono Basin contain several sets of lateral and terminal 
moraines representing multiple stages of glaciation. The semi-arid climate with slow 
weathering rates has preserved sequences of nested younger moraines within older ones. 
There is a well established relative chronology (Sharp and Birman 1963, Burke and 
Birkeland 1979, Gillespie 1982) and recently exposure dating has provide a new set of 
numerical dates (Phillips et al. 1990). The moraines span the late Wisconsin (1 1-25 ka) 
to the Illinoian (130-190 ka) glaciations. We are using the Mono Basin area as a 
"calibration site” to establish remote dating techniques for eventual transfer to the more 
inaccessible but geomorphically and climatically similar moraines of the South American 
Andes Mountains. Planned polarimetric SAR imagery acquired by JPL AIRSAR (South 
American Campaign) and SIR-C (Andes super-site) will be analyzed to establish 
chronologies of previously undated moraine sequences in a study of Pleistocene climatic 
change in the Southern Hemisphere. 

The dry climate and sparse vegetation is also favorable for correlation of ground 
surface roughness with radar polarization signature. The slow weathering processes acting 
over thousands of years reduce the size, frequency and angularity of surface boulders 
while increasing soil development on the moraines. Field observations based on this 
hypothesis (Sharp and Birman 1963) result in relative ages consistent with those inferred 
from nested position within the valley. Younger moraines, therefore, will appear rougher 
than the older smoother moraines at scales measurable at AIRSAR wavelengths. 
Previously documented effects of ground surface roughness on polarization signatures 
(Evans et al. 1988, Ray et al. 1991) suggest that analysis of moraine polarization 
signatures can be useful for relative dating. The technique may be extended to predict 
numerical ages. 

The data set reported here were acquired on September 8, 1989 with the JPL 
AIRSAR collecting polarimetric imagery at C- (5.6 cm), L-(24 cm) and P-band(68 cm) 
with a flight-line parallel to the strike of the mountains. Phase calibration was performed 
on the analyzed scene by setting the co-phase of a smooth lake to zero as described in 
Zebker and Lou (1990). Absolute amplitude calibration was not possible because comer 
reflectors were not deployed. 

2. ANALYSIS 

The Walker Valley moraine complex (Bloody and Saw Mill Canyon) was chosen 
for initial analysis because it contains moraines from five distinguishable glacial episodes 
(Table 1). Estimates of the surface roughness of the crests of the moraines were 
determined from 1-D Fourier transforms of a series of 15 m linear push-rod profiles (10 
cm horizontal spacing, 1.0 mm vertical resolution) from the moraine crests (Fox 1989). 
This measure of surface roughness was found to correlate with the relative moraine ages 
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(Table 1). The crests were measured because they were the areas of minimum soil 
development and least vegetation, enhancing boulder exposure, and thus provide the 
optimal discrimination based on ground surface roughness. Efforts were made to avoid 
vegetation during profiling. Fox found that surface roughness correlated best with the 
image pixel DN values for P-band HH when compared to linear polarizations of the three 
bands (Fox 1991). The following preliminary analysis will focus on P and L-band 
polarization signatures because of their sensitivity to the large scale roughness presented 
by individual boulders. 


Moraine (Glaciation) 
Tioga 
Tenya 

Tahoe Younger 
Tahoe Older 
Mono Basin 


Age [ka] (best estimate) 
14 a 
40 a 
60-85 b 
<118 c 
130 a 


Surface Roughness 
6.70 
6.51 
6.37 
6.22 
5.12 


Table 1. Walker Valley moraine classification with estimated age and measured ground 
surface roughness, a- (Crook and Gillespie 1986) b- (Bursik 1989) c- (Gillespie 1982) 

Polarization signatures were created for each moraine from pixels along and 
adjacent to their crests’. The initial results are based on calculations from typical 
signatures of each moraine. Further research will include more samples for a statistically 
sound study. 

The scattering from an object will become more diffuse as the roughness of the 
target increases relative to the illuminating wavelength. Polarization signatures reflect 
the degree of diffuse scattering in their measure of co- and cross-polarized radar cross 
section as a function of all possible transmit polarizations. The most common way of 
estimating the diffuse component is with the pedestal height (vanZyl et al. 1987) 
calculated as the ratio of the co-polarized minimum to the signature maximum. An 
alternative approach which uses only circular polarization for estimation of diffusivity is 
the circular polarization ratio (referred to here after as CPR) as defined by Ray et al. 
(1991). The right-hand circular polarization ratio is calculated by dividing the cross- 
polarization right-hand response (transmit RH receive LH, expected for a smooth object) 
by the co-polarizied right-hand response (transmit RH receive RH expected from a diffuse 
object), thus comparing the directly reflected power to the diffusely scattered power. 
Rough objects would theoretically have lower CPR’s than smooth objects (relative to 
wavelength). 


Pedestal heights and right-hand CPR (direct/diffuse) values for the five moraines 
at P- and L-bands are shown in Table 2. The P-band data show consistent trends for 
decreasing roughness with successive glaciations. The pedestal height decreases and the 
right-hand CPR increases. The L-band data show the same trends except for the Mono 
Basin moraine, it does not appear to be the smoothest (oldest) using either the pedestal or 
the CPR technique. 


Moraine 

(Glaciation) 

Tioga 

Tenya 

Tahoe Younger 
Tahoe Older 
Mono Basin 



Pedestal Height 

R-H Circular Polarization Ratio 

L 

P 

L 

P 

.1280 

.1633 

.4621 

.1243 

.1230 

.1587 

.5776 

.2143 

.0888 

.0904 

.6549 

.4429 

.0714 

.0702 

1.000 

.6143 

.1178 

.06442 

.7698 

1.000 

Valley 

moraine classification with L- and P-band pedestal heights and 


right-hand circular polarization ratios (normalized) 
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3. SUMMARY 


The preliminary results indicate that the polarization signatures from the glacial 
moraines of the Walker Valley may be used as a relative dating tool. Specifically, the 
parameters of pedestal height and the ratio of circular expected to circular diffuse 
polarization have shown correlation with moraine age at P-band and to a lesser degree at 
L-band. The established correlation of relative moraine age with crest ground surface 
roughness provides the physical link from scattering mechanism to relative age. The 
violation of the trends at L-band for the Mono Basin moraine may be the result of diffuse 
scattering from vegetation creating the “appearance” of a rougher surface. The Mono 
Basin moraines have the most densely vegetated crests. The ground shrubs may have less 
of a response at the longer P-band wavelength. This hypothesis may be further supported 
if the Mono Basin moraine appears even rougher at C-band, 

Additional signatures need to be acquired at more locations representing a diverse 
set of incidence angles, local incidence angles, pixel numbers used for signature, and crest 
orientations to confirm the observed trends. The shape of the polarization signatures will 
be further analyzed for implications of the nature of the scattering mechanism. A more 
detailed and comprehensive characterization of the moraine surface roughness is planned 
through a high resolution DEM to be generated from low altitude stereo photogrammetry. 
The potential exists for relative dating of moraines in inaccessible areas from remote 
polarimertic SAR platforms. 
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1. INTRODUCTION 

The part of the troposphere influenced by the surface of the earth is termed 
the atmospheric boundary layer {Stull 1988). Flow within this layer is influenced by 
the roughness of the surface; rougher surfaces induce more turbulence than smoother 
surfaces and, hence, higher atmospheric transfer rates across the surface. Roughness 
elements also shield crodiblc particles, thus decreasing the transport of windblown 
particles. Therefore, the aerodynamic roughness length (z 0 ) is an important 
parameter in acolian and atmospheric boundary layer processes as it describes the 
aerodynamic properties of the underlying surface, zq * s assumed to be independent 
of wind velocity or height, and dependent only on the surface topography ( Bagnold 
1941; Greeley andlversen 1985). It is determined using in situ measurements of the 
wind speed distribution as a function of height. For dry , unvegetated soils the 
intensity of the radar backscatlcr (a 0 ) is affected primarily by surface roughness at a 
scale comparable with the radar wavelength. Thus, both wind and radar respond to 
surface roughness variations on a scale of a few meters or less. Greeley cl al. {1991) 
showed the existence of a correlation between zq and cr®. This correlation was based 
on measurements over lava flows, alluvial fans, and playas in the southwest deserts 
of the United States. In this report we show that the two parameters behave similarly 
also when there are small changes over a relatively homogeneous surface. 

2. EXPERIMENTAL DETAILS 

A study was conducted during the Spring of 1991, in which the atmospheric 
boundary layer was assessed and compared with a 0 values obtained from AIRSAR 
data. Five masts equipped with wind anemometers, wind vanes, and temperature 
sensors were deployed over an unvegetated alluvial fan on the eastern edge of the 
Stove Pipe Wells Dune Field in Death Valley, California. The alluvial fan lies al sea 
level below the Grapevine Mountains, and receives its deposits primarily from Mud 
Canyon. The predominant airflow in this area is NNW-SSE, across the fan. 

Micromctcorologic data were collected for a period of 10 weeks and 
reduced following a procedure by While ( unpublished) to derive the surface shear 
velocity, the Richardson number {Stull 1988), z 0 per wind profile, and a mean z 0 per 
data set. Each data set consisted of approximately 100 wind profiles for a given wind 
direction per site (Table 1). AIRSAR data were obtained at the site using bands C 
(5.6 cm), L (24 cm), and P (68 cm). The radar data were calibrated using POLCAL 
software {vanZyl el al. 1990), based on trihedral comer reflectors that were deployed 
prior to the imaging. 
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3. RESULTS 


Mean o° values were extracted for the location of each of the five 
masts (Table 2) (values did not need to be normalized for variations in incidence 
angles, as all incidence angles were between 45° to 47°). Figure 1 shows that 
and zq change simultaneously and in the same direction for winds from the 
north-northwest, suggesting that they are responding to similar surficial changes. 
These variations are probably induced by roughness changes comparable with 
the radar wavelength. Surface roughness measurements were conducted using 
an electronic distance meter, a surface template, and a laser profiler. The surface 
roughness data are being analyzed to determine if the changes in surface 
roughness corroborate the radar and wind results. Winds from the south- 
southeast showed larger variations of z 0 , suggesting that they are affected by the 

fluvial dissection found south of the study site. 

4. CONCLUSIONS 

Based on the north-northwesterly winds we conclude that Zq and a 0 

have similar spatial patterns even when measured over a relatively homogeneous 
surface. L-band cross-polarized data correlate best with z 0 over a dry, 

un vegetated alluvial fan. Further work is required to determine the scales of 
roughness affecting these parameters. 


Table 7. zq across the alluvial fan 


Mast Distance from 

northwest to 
southeast (m) 
2 200 

3 400 

4 700 

5 1000 


7-o (m) 

southeasterly 

winds 

0.00608 

0.00055 

0.00806 

0.00240 


z 0 (m) 

northwesterly 

winds 

0.00208 

0.00076 

0.00288 

0.00310 


Table 2. cP variations across the alluvial fan 


Mast C-band (dB) L-band (dB) P-band (dB) 



HH 

HV 

VV 

HH 

HV 

VV 

HH 

HV 

VV 

1 

-12.72 

-2,0.97 

-12.32 

-21.92 

-34.70 

-20.06 

-31.31 

47.01 

-27.14 

2 

-14.29 

-22.40 

-12.77 

-24.43 

-36.94 

-22.46 

-34.27 

49.73 

-30.78 

3 

-13.87 

-22.52 

-12.85 

-25.34 

-40.61 

-22.96 

-35.00 

-50.65 

-31.87 

4 

-13.16 

-22.47 

-13.00 

-23.40 

-37.39 

-21.18 

-32.74 

48.78 

-29.67 

5 

-12.35 

-21.71 

-11.66 

-22.36 

-35.66 

-20.48 

-32.20 

49.01 

-29.04 


44 







o° L-HH ~°~ z a 



Figure 1 : Mean radar backscatlcr in L-band (24 cm) , and zq values for northwesterly 
winds versus the locations of the masts across the alluvial fan. For the most 
northern station no z 0 value was determined due to instrument problems. 
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L INTRODUCTION 

Polarimetric Airborne Synthetic Aperture Radar (AIRSAR) data were collected 
for the Geologic Remote Sensing Field Experiment (GRSFE) over Death Valley, 
California, USA, in September 1989 (Evans and Arvidson, 1990; Arvidson et i, 1991). 
AIRSAR is a four-look, quad-polarization, three frequency instrument. It collects 
measurements at C-band (5.66 cm), L-band (23.98 cm), and P-band (68.13 cm), and has a 
GIFOV of 10 meters and a swath width of 12 kilometers. Because the radar measures at 
three wavelengths, different scales of surface roughness are measured. Also, dielectric 
constants can be calculated from the data (Zebker et al, 1987). 

The AIRSAR data were calibrated using in-scene trihedral comer reflectors to 
remove cross-talk; and to calibrate the phase, amplitude, and co-channel gain imbalance 
(van Zyl, 1990). The calibration allows for the extraction of accurate values of rms 
surface roughness, dielectric constants, a 0 backscatter and polarization information. The 
radar data sets allow quantitative characterization of small scale surface structure of 
geologic units, providing information about the physical and chemical processes that 
control the surface morphology. Combining the quantitative information extracted from 
the radar data with other remotely sensed data sets allows discrimination, identification and 
mapping of geologic units that may be difficult to discern using conventional techniques. 

2. INVERSION AND ANALYSIS 

The first-order small perturbation model (van Zyl ct al ,1991; Valenzuela, 1967) 
was used to estimate the surface power spectral density at every pixel by performing an 
inversion of the AIRSAR data. This model is valid only for very smooth surfaces. 
Therefore, only playa, smooth salt pans and alluvial fan surfaces were used in this study. 
Results from the small perturbation model inversion are three values, one for each of the 
radar frequencies, that describe the power spectral density of the surface. The power 
spectrum of a geologic surface is approximately linear in log-log space (Brown and 
Scholz, 1985). Fitting the three points from the inversion with a line using a least- 
squares method produces slope and intercept values that allow calculation of the fractal 
dimension of the surface and a rms surface roughness value. 

The slope of the power spectrum is related to the two-dimensional fractal 
dimension of the surface. The fractal dimension of a surface describes the scaling 
properties of the topography (Mandelbrot, 1982). A surface may have a fractal dimension 
between 2 and 3 and as the fractal dimension increases, heights of nearby points become 
more independent (Brown and Scholz, 1985). The intercept of the power spectrum can be 
directly related to a rms surface roughness using forward modelling. Using the fractal 
dimension and rms surface roughness calculated from the radar inversion power spectrum, 
a synthetic three dimensional plot can be made that represents the surface (Huang and 
Turcotte, 1989; Kierein-Young etal, 1992). 
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3. RESULTS AND CONCLUSIONS 


Radar inversions were performed using two dielectric constant cases for polygon 
averages of four surfaces in Death Valley. The surfaces include a playa, smooth salt, 
alluvia] fan with desert pavement, and a rough alluvial fan. In both of the dielectric 
constant cases, the rms surface roughness is smallest for the playa, increases for the 
smooth salt and fan with desert pavement, and is largest for the rough fan. The fractal 
dimension is largest for the smoothest surface, the playa, and smaller for the rougher 
surfaces. This indicates that nearby points on the smoother surfaces are more independent 
and nearby points on the rougher surfaces are more correlated. Therefore, the smooth 
surfaces have a more equal distribution of all wavelengths of roughness while the rough 
surfaces have more longer wavelength roughness. Figure 1 shows the synthetic surfaces 
for the four areas with the dielectric constants equal to 3.0. In the case when dielectric 
constants are calculated in the inversion, some bad values were obtained so more analysis 
is needed. In this case the fractal dimensions are slightly smaller and the rms roughness 
values are higher than in the dielectric constant equals 3.0 case. The values of the 
dielectric constants calculated in the inversion fall within a reasonable range for the 
surfaces studied (Ulaby et al, 1982). Dielectric field and laboratory measurements will be 
made to verify the accuracy of these values. 

Combining knowledge of the surface structure and roughness of geologic units 
with data from other instruments allows for the determination of the processes that 
control the surface. For example, water vapor and liquid water maps were generated from 
Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data in Death Valley (Kierein- 
Young and Kruse, 1991). These maps show the distribution of a hydrological process, 
the water content of the surface. The water distribution was found to correspond with the 
surface roughness variations where wet areas are smooth and drier areas rougher. The 
combination of surface water distribution with surface structure information will help in 
determining why the salt flats are spatially variant. 
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Figure 1. Synthetic surfaces calculated from radar inversion results assuming dielectric constants equal to 3.0. 
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Abstract 

We have augmented the NASA DC-8 AIRSAR instrument with a pair of C-band antennas displaced across 
track to form an interferometer sensitive to topographic variations of the Earth’s surface. The antennas were 
developed by the Italian consortium Co.Ri.S.T. A., under contract to the Italian Space Agency (ASI), while 
the AIRSAR instrument and modifications to it supporting TOPSAR were sponsored by NASA. A new 
data processor was developed at JPL for producing the topographic maps, and a second processor was 
developed at Co.Ri.S.T. A. All the results presented below were processed at JPL. During the 1991 DC-8 
flight campaign, data were acquired over several sites in the United Slates and Europe, and topographic 
maps were produced from several of l fuse flight, lines. Analysis of the results indicate that statistical errors 
are in the 2-8 m range for flat terrain and in the 4-5 m range for mountainous areas. 

Introduction 

Wc have implemented an interferometric synthetic aperture radar system for topographic mapping ap- 
plications (TOPSAR) on the NASA DC-8 aircraft. NASA /JPI, currently operates a multi frequency (P, 
L, and C bands), multipolarimctric radar (AIRSAR) on board this aircraft. The TOPSAR implementa- 
tion uses much of the existing AIRSAR hardware, although several modifications were required to optimize 
performance in the topographic mapping mode. When in use, TOPSAR effectively replaces the C-band 
polarimeter instrument, but the remaining L- and P-band systems are undisturbed and operate together 
with the topographic mapper* therefore the combined instrument produces simultaneous L- and P-band 
hilly polarirnetric plus C-band VV polarization backscaller images in addition to the topographic product. 

In this paper we first describe our implementation, including system design parameters, required mod- 
ifications to the existing AIRSAR hardware, and new antenna design. We then discuss data processing 
strategy and display some example interferometric images. Finally, we analyze an image acquired over Ft. 
Irwin, California, in terms of height accuracy. 

Background 

Interferometric radar has been been proposed and successfully demonstrated as a topographic mapping 
technique by Graham [1], Zebker and Goldstein [2], and Gabriel and Goldstein [3]. Two distinct imple- 
mentation approaches have been discussed Tor topographic radar interferometers which differ in how the 
interferometric baseline is formed. In one case a single antenna and radar system illuminates a given surface 
at two different times but with nearly the same viewing geometry, forming a synthetic interferometer. This 
case has been implemented both from spaceborne platforms (see, for example, Goldstein el al., [4]) and also 
from aircraft (Gray and Farris- Manning [5]). In the second approach the baseline is formed by two physical 
antennas which illuminate a given area on the ground simultaneously. I his is the approach used originally 
by Graham [1] and also by Zebker and Goldstein [2] for the NASA CV-990 radar; it is the approach used 
here. 

Review of theory 

Consider a set of two antennas Al and A2 as shown in figure 1. T’he surface topography is given by z(y), 
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h is the aircraft altitude, the baseline distance is B, the range to a point on the ground is /?, the look angle 
9, and the angle of the baseline with respect to horizontal is a. Radar echoes are transmitted from antenna 
A1 and received simultaneously at Al and A2, thus the difference in path lengths is 6, which depends on the 
baseline distance, baseline angle, range, look angle, and the height of the point z(y). The measured phase of 
the interferometer is directly proportional to this distance, with the constant of proportionality . A little 
algebra and geometry yield the following equations for height as a function of these parameters: 

- e) - ( " +<> t /~ g2 p> 

z(y) ~ h - p cos a cos(a - 9) - p sin a sin(a - 9) (3) 

where 4> is the measured phase, and A is the wavelength. 

Differentiation of (1-3) yields the error in height estimate as a function of the error in phase estimate: 

° h = 27 rT? a “ cos a tan ( a ” 0)] a 4> (4) 

where and are the standard deviations of height and phase, respectively. 

There exists an optimum baseline distance that minimizes the total height error [Rodriguez and Martin, 
6]; this optimum is quite broad and for high (> 10 dB) signal to noise ratio systems can be anywhere in the 
range of 0.2 - 0.8 of the critical value which may be calculated by ( Zebker and Villasenor [7]): 

B c = tan 9 (5) 

where Ry is the slant range resolution. Note that this equation differs from that of Zebker and Villasenor by 
a factor of 2 cos# sin# as they were restricted to horizontal baselines, used the ground rather than slant range 
resolution, and assumed a single antenna in repeat-track configuration rather than two physical antennas as 
we use in TOPSAR. Another significant error source results from errors in knowledge of the aircraft attitude, 
the most important component of which is roll angle, which in our system translates directly to an error in 
look angle. Again differentiation yields 

= #sin 9a$ (6) 

Thus, good motion compensation is a prerequisite of any practical system. 

TOPSAR design 

Several TOPSAR system parameters are summarized in Table 1. 

Table 1. TOPSAR radar system parameters 


Parameter 


TOPSAR value 


Peak power, watts 

1000 

Pulse rate, Hz 

600 nominal 

Pulse length, fi sec 

5.0 

Antenna length, m 

1.6 

Antenna width, m 

0.11 

Antenna gain, dB 

25 

Range bandwidth, MHz 

40 

Receiver noise temperature, K 

2100 

Antenna baseline, m 

2.58 

Baseline angle (a), deg 

62.77 


Given the baseline separation and the signal to noise ratio we can then analyze performance of the 
interferometer. Figure 2 ( Li and Goldstein [8], Zebker and Villasenor [7]) gives the expected phase error as 
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a function of signal to noise ratio and number of looks for a radar interferometer. Thus we may expect a 
phase error of 3.3° if the data are processed to 16 looks assuming an SNR of 13 dB, which implies a statistical 
error of 1.50 m in height precision. Knowledge of the roll angle is also needed for accurate data reduction. 
At present our system can correct for attitude errors to the 0.015° level, which (6) implies a height error of 
1.9 m. 

Figure 3 shows a set of data acquired over Fort Irwin in the Mojave Desert in California. The data shown 
represent ground coordinate rectified topographic maps in which the color contour interval is 6 m, thus one 
trip around the color wheel corresponds to 96 m for our 16 entry color table. The spatial resolution is 10 m. 
The accuracy of the map obtained over Ft. Irwin is described in the next section. 

Verification procedure 

Verification of this mapping technique involves comparing the radar-derived height map with pre-existing 
high- resolution digital elevation models (DEMs) and characterizing the differences. Comparison was done 
on a data set from the Ft. Irwin area in California using the USGS 7.5’ DEM with 30 m grid spacing and 
7 m rms height error. Figures 4 shows a sample cut through the DEM and radar-generated terrain maps, 
and the difference between them is shown in the upper part of the figure. It is clear from examination of 
this figure that the error is dominated by differences in areas of high relief, and is caused both by errors in 
overall cross-track slope and by position errors. 
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1. TOPSAR geometry. Radar echoes are transmitted from an- 
tenna A1 and received simultaneously at Al and A 2. The phase 
difference of the two echoes is proportional to the difference in 
path lengths S t which depends on the baseline distance B, base- 
line angle a, range p, look angle 9, aircraft altitude h, and the 
height of the point z(y). 



2. Sensitivity of phase standard deviation to signal- to- noise ratio 
and number of looks in processor. Increasing number of looks is 
an effective means to reduce statistical variation, especially for 
the first 8 looks or so. 
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Along-track distance (m) 


3. Topographic data acquired over Fort Irwin in the Mojave 4. DEM / radar map height difference from the Fort Irwin data. 

Desert in California. The data shown represent ground coor- 
dinate rectified topographic maps in which the color contour 
interval is 6 m. 
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This paper presents the preliminary analysis of the C-band cross- 
track interferometric data (XT1) acquired during the MAC Europe 1991 
campaign over the Matera lest site, in Southern Italy. 23 passive 
calibrators (Corner Reflector, CR) and 3 active calibrators (Active Radar 
Calibrator, ARC) were deployed over an area characterized by homogeneous 
background (Fig. 1). Contemporaneously to the flight, a ground truth data 
collection campaign was carried out {COMS.T.A, 1991). 

The DC- 8 NASA airplane was equipped with the JPL A1RSAR, which 
had been modified at the Ames Research Center to accommodate two C- 
band antennas designed and developed by the C0.R1.S.T.A. Consortium and 
to receive the two channels by using the existing hardware (TOPSAR, Zebker 
etaL 1992). 

The research activity has been focused on the development of 
motion compensation algorithms, in order to improve the height 
measurement accuracy of the TOPSAR system. Moreover, we intend to 
derive the system impulse response by comparing and integrating a 
theoretical model, the available ancillary data, and the CRs phase history. 

In the following we present a status report of our research, which 
has been mainly devoted to compute a reference function for azimuth 
compression, that accounts for phase errors due to the airplane dynamics. 
In this work, after the characterization of the doppler histories of the CRs 
and the ARCs, we carry out a correlation between the attitude data and the 
phase of the range-compressed calibrators. Experimental results show a 
satisfactory correspondence between the airplane dynamics described in 
the ancillary data and the one reconstructed from the phase history of the 
calibrators (shown in Fig 2). At this stage we are able to derive a range- 
dependent azimuth reference function which compensates also velocity and 
altitude perturbations. In order to validate this procedure, we perform an 
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image quality analysis for several calibrators, evaluating their one- 
dimensional 3dB resolutions, ISLR (Integrated Sidelobe Ratio), PSLR (Peak 
Sidelobe Ratio) and SCR (Signal — to- Clutter Ratio) [JPL S/R-C Team et aL 
1990), and consequently an estimate of their radar cross section. A 
preliminary comparison of numerical results with theoretical performance, 
ground truth and anechoic chamber data is encouraging, and allows a 
refinement of the azimuth filter. 

Then, the raw data gathered by the two antennas are processed 
separately and combined to form an inlerferogram. The reconstructed 
airplane dynamics is also used to derive the baseline time variation, that is 
necessary for the 2tt phase unwrapping procedure. To this end we have 
developed a computer code based on a two-dimensional algorithm which 
uses the calibrators as ground control points, both to solve the first 2w 
ambiguity and to overcome decorrelation problems when possible. 

Finally, the computed terrain elevation is compared with the 
1:25000 scale contour levels of the Italian Islituto Geografico Militare (1GM), 
in order to assess the height accuracy of the system and to evaluate the 
improvement consequent to a more sophisticated interferometric 
processing. 

This work has been carried out under the sponsorship of the 
Italian Space Agency (ASI). 
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Fig. 1. Location of CRs and ARCs in the test site. 



AZIMUTH SAMPLES #10 


Fig. 2. Differences of the airplane-CR distance along two adjacent range bins in the 

range-compressed data. 




N94-16751 » 

: j 

FITTING A THREE-COMPONENT SCATTERING 
MODEL TO POLARIMETRIC SAR DATA 


A. Freeman and S. Durden 


Jet Propulsion Laboratory 
M/S 300-235, 4800 Oak Grove Drive 
Pasadena, CA 91 109 


1.0 INTRODUCTION 

A major problem in analyzing polarimetric SAR data such 
as that produced by the NASA/JPL AIRSAR is in understanding 
the scattering mechanisms which give rise to features in the 
different polarization parameters. Researchers, on examining 
some polarimetric SAR data from their scene of interest for the 
first time, often notice unusual bright or dark features when 
displaying one of the many possible polarization representations 
of the data (e g. total power, HH, VV or HV cross-section, HH-VV 
phase difference, HH-VV correlation coeffficient, etc.) It is 
usually pertinent to ask what scattering mechanisms give rise to 
the unusual features? 

Much excellent work has recently been done on modelling 
polarimetric radar backscatter for both naturally occurring 
terrain and man-made objects. These models are usually complex, 
and require a large number of input parameters to successfully 
predict the observed backscatter. For example, in modelling the 
backscatter return from a forest, measurements of tree heights 
and diameters, tree density, leaf size and angular distribution, 
branch size and angular distribution, trunk dielectric constants, 
ground roughness and dielectric constant are commonly 
required as inputs. All these models solve the forward problem', 
in predicting backscatter from a number of ground-based 
measurements of the imaged objects. It is difficult, if not 
impossible, to invert these models to provide a unique solution, 
simply because the number of input parameters (the ground 
truth ) is often much larger than the number of output 
parameters (the radar measurements) in the forward problem. 

In this paper, a new technique for fitting a three- 
component scattering mechanism model to the polarimetric SAR 
data itself, without utilising any ground truth measurements, is 
presented. The three scattering mechanism components included 
in the model are volume scatter from randomly oriented dipoles, 
first-order Bragg surface scatter and a dihedral scattering 
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mechanism for two surfaces with different dielectric constants. 

The model fit yields an estimate of the contribution to the total 
backscatter of each of the three components. The backscatter 
contributions can also be compared to give the relative 
percentage weight of each. The model fit has an equal number of 
input parameters (the polarimetric radar backscatter 
measurements) and output parameters (the backscatter 
contributions from each of the three components and two 
parameters describing them). The model can be applied to entire 
images or to small areas within an image to give a first-order 
estimate of the relevant scattering mechanisms. The model has 
been applied to many C-, L- and P-band AIRSAR images of 
different types of terrain. Results will be presented at the 
workshop. 

2. THE MODEL 

The model fit includes three very simple scattering 
mechanisms. First, for volume scattering, it is assumed that the 
radar return is from randomly oriented, very thin cylinder-like 
scatterers. By making several simplifying assumptions the 
second-order statistics of the resulting scattering matrix can be 
derived. After normalizing with respect to the VV cross section, 
these are: 

|shh| 2 ) = |svv| 2 } = 1 , (ShiiSvv) - |shv| 2 ) =1/3 ( 1 ) 

The double-bounce scattering component is modelled by 
scattering from a dihedral corner reflector, where the reflector 
surfaces can be made of different dielectric materials, 
corresponding to a ground-trunk interaction for forests, for 
example. The vertical surface (e g. the trunk) has Fresnel 
reflection coefficients Rth and Rtv for horizontal and vertical 
polarizations, respectively. The horizontal surface (the ground) 
has Fresnel reflection coefficients Rgh and Rgv • The scattering 
matrix for double-bounce scattering is then: 

_ -Rgv Rtv 0 ( 2 ) 

0 Rgh Rth 

The second-order statistics for double-bounce scattering, after 
normalization with respect to the VV term, are: 

|shh| 2 } = N 2 , fsvv| 2 ) = 1 , (shhS vv) = -a , |shv| 2 ) = 0 

where a = Rgh Rth / Rgv Rtv (3) 

For the surface scatter, a first-order Bragg model is used, with 
second order statistics (after normalization): 

|shh| 2 ) = |(3| 2 . \ s vv| 2 ) - 1 . (shhSvv) = P , |shv| 2 ) = 0 (4) 
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where (3 is real For all of these backscatter components, it is 
assumed that the like- and cross-polarized returns are 
uncorrelated, and that the backscatter is reciprocal. Further, 
assuming that the volume, double-bounce and surface scatter 
components are uncorrelated, the total second order statistics are 
the sum of the above statistics for the individual mechanisms. 
Thus the model for the total backscatter is: 

(ShhI 2 ) - A Iff * /-^IPl 2 + fv 


|sw| 2 ) = fs + fa * fv 


(ShhSvv) = /* P - fd OL + fy / 3 


js H v| 2 ) = fy/ 3 


(5) 


where f s , fd and f v are the surface, double-bounce and volume 
scatter contributions to the VV cross section. 

The above model gives four equations in five unknowns. A 
solution can only be found if one of the unknowns is fixed. Since 
neither the surface nor the double-bounce mechanism 
contribute to the HV term, this can be used to estimate the volume 
scatter contribution directly. This can then be subtracted off the 
three remaining terms, leaving three equations in four 

unknowns. The sign of the real part of the residual ShhSvv term 
can then be used to decide whether double-bounce or surface 

scatter is dominant. If Re (ShhSvv) is positive, then surface scatter 

is determined to be dominant, and (X is set to 1. If Re (ShhSvv) is 
negative, then double-bounce scatter is determined to be 
dominant, and (5 is set to 1. Then fs , fd • fv and P or a can be 
estimated from the residual radar measurements. 

3. SUMMARY 

A new technique for fitting simple backscatter 
mechanisms to polarimetric SAR data has been presented. The 
model can be used to determine to first order what the dominant 
scattering mechanisms are which give rise to observed 
backscatter in polarimetric SAR data. Results of application of the 
model to AIRSAR images of different types of terrain will be 
presented at the workshop. 
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1. INTRODUCTION 

Based on symmetry properties, polarimetric remote sensing of geo- 
physical media is studied in this paper. From the viewpoint of symmetry 
groups, media with reflection, rotation, azimuthal, and centrical symme- 
tries are considered. The symmetries impose relations among polarimet- 
ric scattering coefficients, which are valid to all scattering mechanisms 
in the symmetrical configurations. Various orientation distributions of 
non-spherical scatterers can be identified from the scattering coefficients 
by a comparison with the symmetry calculations. Experimental observa- 
tions are then analyzed for many geophysical scenes acquired with the Jet 
Propulsion Laboratory (JPL) airborne polarimetric SAR at microwave 
frequencies over sea ice and vegetation. Polarimetric characteristics of dif- 
ferent ice types are compared wuth symmetry behaviors. The polarimetric 
response of a tropical rain forest reveals characteristics close to the cen- 
trical symmetry properties, which can be used as a distributed target to 
relatively calibrate polarimetric radars without any deployment of man- 
made calibration targets. 

2. SYMMETRY PROPERTIES 

Relations among polarimetric backseat tering coefficients have been 
derived (Nghiem et al. 1992) for both reciprocal and non-reciprocal geo- 
physical media w r ith reflection, rotation, azimuthal, and centrical symme- 
tries. The corresponding symmetry groups can be constructed from three 
fundamental operations: mirror reflection, axja) rotation, and linear trans- 
lation (H amermt sh , 1972). The derivations are based on the invariance 
of the scattering coefficients under the symmetry transformations and ap- 
plicable to all scattering mechanisms in the symmetrica! configurations. 
These symmetries impose on the scattering coefficients a number of equa- 
tions. which reduce the number of independent parameters in the covari- 
ance or Mueller matrix as summarized in Table 1. Centrical symmetry 
inherits azimuthal symmetry characteristics at all incident angles. 


Table 1: Number of Independent Parameters in Covariance 
Matrix or Mueller Matrix 



Non-reciprocal 1 

Reciprocal 

No symmetry 

16 

9 

Reflection Symmetry 

8 

- 

5 

Rotation Symmetry 

6 i 

- • -t 

3 

Azimuthal Symmetry 

i — — — - — 

3 • 

. . - , -1 - 

~2 
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3. ORIENTATION DISTRIBUTIONS 


Non-spheriral scatterers can take on needle-like shape or disc- 
like shape. The distribution of scatterer orientations is dependent on the 
scatterer species, environmental effects, and physiologic conditions. For 
species with no azimuthal preference such as leaves in trees, several distri- 
butions have been reported including spherical, uniform, planophile. pla- 
giothile, erectophile, or extremophile distributions (Gotl and Stnbd 198J). 
Probability density functions of orientations are given in terms of inclina- 
tion angles 6 from 0 to 90°. These orientation distributions have certain 
symmetries which affect polarimetric scattering coefficients. In each case 
of the distributions, the ratio e — a hv )a hh is calculated under independent 
scattering assumption from a theoretical model (Tsang et al. 1985) and 
then compared with the quantity eo = (1 — \p )/2 as suggested by the cen- 
trical symmetry where p is the correlation coefficient, between the hh and 
vv returns. For centrical symmetry, e = at arbitrary incident angles. 
Referred to this symmetry, the deviation of e from cq therefore indicates 
how the scatterers are structured. The results for the aforementioned ori- 
entation distributions show that, in general, c/c 0 is larger for more ver- 
tical orientation regardless of scatterer shape. For spherical distribution, 
e / e o - ] manifests the centrical symmetry. 

4, APPLICATIONS TO JPL SAR DATA 

In this section, an analysis based on the symmetry is applied to 
JPL SAR data for sea ice and vegetation. Polarimetric SAR scenes of 
sea ice in the Beaufort sea were acquired during March 1988 in a series of 
NASA DC-8 airborne flights (Cavalieri, 1988). Weather and sea ice char- 
acteristic data were collected continuously during March at the Applied 
Physics Laboratory's drifting ice station (APLJS’88) located north-east 
of Prudhoe Bay, Alaska. Meteorological reports and ire measurements 
were recorded (Wen d al 1989). Sea ice conditions consist of a transition 
from the extensive coastal first-year ice region to multi-year ice pack to 
the north. Air temperatures varied between 12" and 18 on M arcli 1 1 
and colder (below -25 c ) during an earlier period in March. First-year ice 
in the vicinity of the ice station was 1.5-m to 2.4-in thick and covered by 
a snow layer of 15.0-cm average depth. Multi-year ice was also covered by 
snow and hummocked up to 6.0-m height. Ice divergence occurring close 
to the period of dat a collection was experienced on March 10 when high 
ice drift velocities were recorded. Sea ice scenes on March 11 reveal im- 
ages of new cracks and leads created by the divergent ice motion. Leads 
rapidly frozen under the cold condition became new thin ice formations. 
The results indicate that first-year ice signatures quite depart from centri- 
cal symmetry behavior as compared to multi-year ice. Structurally, first- 
year ice has a preferential vertical orientation of brine inclusions (Weeks 
and Ackley 1982) while scatterers in multi-year ice are more randomly dis- 
tributed. Comparisons of the measured values for the polarimetric returns 
from sea ice with the symmetry calculations reveal the structural informa- 
tion which helps identify the ice types. 

For vegetation. SAR scenes acquired over various forest types, in- 
cluding pine, mixed, and tropical rain forests, are considered. In the 1989 
MAESTRO- 1 Campaign, an experiment was conducted in Les Landes 
pine forest, south-west France, where the JPL SAR data were collected to 
relate SAR measurements and forest biophysical parameters (he Toan et 
al. 1991). This scene consists of maritime pine trees at different ages from 
young to mature stages (more than 40 years) Also in June 1989, SAR 
images were obtained for Mt. Shasta in northern California Forest ar- 
eas in Mt. Shasta are dominated by ponderosa pine and white fir species 
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( Ztbker it al. 1991). In Belize, a tropical rain forest was imaged. The for- 
est in this scene has dense canopy and understory with various tropical 
species. At C-band frequency, the relation e - f () . small phase of p, and 
the balance between hh and vti returns are well satisfied at all incident 
angle, especially for the Belize forest. This can be explained based on the 
symmetry properties. Azimuthal symmetry is often observed on leaves 
orientation in forest canopies. For elevation orientations of vegetation el- 
ements, the distributions have been discussed in the last section. Foliage 
composed of leaves with spherical orientation distribution has centrical 
symmetry wdiich requires the polanmetric scattering coefficients to satisfy 
e £(1 , Ini p -- 0, and 7 - a vv ja hh 1 at arbitrary incident angles. These 
relations are the direct consequence of the* centrical symmetry due to the 
random orientation of leaves (spherical distribution). When the frequency' 
is low such that the electromagnetic wave can penetrate through the fo- 
liage canopy, the centrical symmetry can be destroyed due to the horizon- 
tal branches or other structures of tree elements. From natural distributed 
largets, reflection together with centrical symmetries will allow the full 
relative calibration, including the cross talk and the channel imbalance, to 
all incident angles without tfie deployment of man-made calibration tar- 
gets. 
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1. INTRODUCTION 

A new technique is presented in this paper for calibrating polarimetric 
synthetic aperature radar (SAR) images using only the responses from natural 
distributed targets. The model for polarimetric radars is assumed to be (Yueh et 
ah, 1990) 

X = cRST (1) 

where X is the measured scattering matrix corresponding to the target scatter- 
ing matrix S distorted by the system matrices T and R (in general T ^ R l ). 

To allow for the polarimetric calibration using only distributed targets and cor- 
ner reflectors, van Zyl (1990) assumed a reciprocal polarimetric radar model with 
T - R*] when applied for JPL SAR data, a heuristic symmetrization pro- 
cedure is used by POLCAL to compensate the phase difference between the mea- 
sured HV and VH responses and then take the average of both. This heuristic 
approach causes some non-removable cross-polarization responses for corner re- 
flectors, which can be avoided by a rigorous symmetrization method based on 
reciprocity (Yueh et al., 1991). After the radar is made reciprocal, a new algo- 
rithm based on the responses from distributed targets with reflection symmetry 
is developed to estimate the cross-talk parameters (Yueh et al., 1992). The new 
algorithm never experiences problems in convergence and is also found to converge 
faster than the existing routines implemented for POLCAL. When the new tech- 
nique is implemented for the JPL polarimetric data, symmetrization and cross- 
talk removal are performed on a line-by-line (azimuth) basis. After the cross-talks 
are removed from the entire image, phase and amplitude calibrations are carried 
out by selecting distributed targets either with azimuthal symmetry along the 
looking direction or with some well-known volume and surface scattering mech- 
anisms to estimate the relative phases and amplitude responses of the horizontal 
and vertical channels. 

2 APPLICATIONS TO SAR DATA 

After this new technique is applied to the L-band polarimetric data col- 
lected over the Mt. Shasta area in 1989, the responses of corner reflectors were 
extracted after each step to assess the effectiveness of the methods. Tables 1 and 
2 show the covariance matrix parameters of the corner reflectors after cross-talk 
removal. The data presented in Table 1 are symmetrized by the POLCAL’s ap- 
proach before cross-talks are removed, whereas the data shown in Table 2 are 
symmetrized based on reciprocity. Comparing Tables 1 and 2 for both e and 7 
indicates that the new technique yields superior performance. Note that the same 
algorithm for cross-talk removal is applied for both cases; therefore the difference 


62 


is caused by the symmetrization procedure alone. 


Trihedral 


7 

P 

0 


CRl 

-16.78 

-1.34 

0.997Z31.6 

0.92Z145 

0.95Z- 113 

CR2 

-21.46 

-0.97 

0.995Z29.9 

1.04Z144 

1.05Z- 115 

CR3 

-24.79 

-1.67 

0.991Z31.9 

0.96Z114 

0.98Z - 86 

CR4 

-16.04 

-1.2 

0.958Z35.5 

0.97Z133 

1.02Z-96 

CR5 

-20.27 

-2.41 

0.983Z30.3 

1.02Z127 

1.04Z - 99 

CR6 

-13.59 

-3.01 

0.985Z22.1 

1.01 Z 126 

1.00Z- 101 

CR7 

-13.44 

-1.68 

0.99Z22.1 

0.99Z124 

1.02Z- 102 

CR8 

-19.64 

2.40 

1 .00Z 10.9 

0.85Z — 31.5 

0.86Z40.6 


Table 1. Covariance parameters of trihedral reflectors after symmetrization (POL- 
CAL) and cross- talk removal, e = \HV\/\H H\, y — ^yV\/\HH\, p\ corre- 
lation of HH and VV, 0: correlation of HH and HV , and £: correlation 
of HV and VV. e and y are in dB. The phases of p, 0, and £ are in de- 
grees. 


Trihedral 

e 

7 

P 

0 

£ 

CRl 

-24.0 

0.07 

0.994Z80.6 

0.76Z176 

0.79Z - 100 

CR2 

-100 

0.00 

0.991Z81.6 

- 

- 

CR3 

-100 

-0.96 

0.984Z83.8 

- 

- 

CR4 

-24.07 

-0.07 

0.969Z84.6 

1.01Z137 

1.15Z-42 

CR5 

-100 

-1.38 

0.992Z81.2 

- 


CR6 

-100 

-0.34 

0.992Z69.6 

- 


CR7 

-100 

0.21 

0.998Z68.7 

- 

- 

CR8 

-100 

0.48 

0.991Z64.0 

- 

- 


Table 2. Covariance parameters of trihedral reflectors after symmetrization and 
cross-talk removal based on the new technique. 


For the phase calibration, we determined the relative phase errors be- 
tween H and V channels for nine arbitrarily selected areas in the image, and lin- 
early fit the phase correction factors along the range. Each azimuthal line is then 
calibrated by the interpolated phase factor. The resulting responses of corner re- 
flectors are shown in Table 3. Comparing the phases of p between Tables 2 and 
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3, we can see that the relative phase errors between HH and VV are greatly re- 
duced. 


From the 7 of the corner reflectors, it seems that the amplitudes of H 
and V channels are well balanced for L-band. Thus, additional correction on 
channel imbalance does not seem to be necessary. 
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Trihedral 

e 

7 

P 

P 

£ 

CRl 

-24.0 

0.07 

■ama 

0.76Z131 

0.79Z — 145 

CR2 

-100 

0.00 

1 

- 

- 



-0.96 

0.984Z — 6.5 

- 

- 

CR4 

-24.07 

-0.07 

0.969Z-3.3 

1.01Z94 

1.15Z-86 

CR5 

-100 

-1.38 

0.992Z- 7.0 

- 

- 

CR 6 

Ha 

-0.34 

0.992Z - 10.5 

- 

- 

CR7 


0.21 

0.998Z - 11.8 

- 

- 

CR 8 

-100 

0.48 

0.991Z- 16.7 

- 

- 


Table 3. Covariance parameters of trihedral reflectors after symmetrization, cross- 
talk removal, and phase calibration based on the new technique. 
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1. INTRODUCTION 

Microwave remote sensing technique has a high 
potential for measuring soil moisture due to the 
large contrast in dielectric constant of dry and 
wet soils (Ulaby et al. 1982). Recent work by 
Pults et al. (1990) demonstrated the use of X/C- 
band data for quantitative surface soil moisture 
extraction from airborne SAR system. Similar 
technique was adopted in this study using 
polarimetric SAR data acquired with the JPL-AIRSAR 
system over the Mahantango watershed area in 
central Pennsylvania during July 1990. The data 
sets reported here include C-, L- and P-bands of 
July 10, 13, 15 and 17, 1990. 


2. RESULTS AND DISCUSSION 

SAR data was calibrated and registered with 
digital elevation data using a second degree 
polynomial to study the possible effects of local 
topography on averaged over a corn field (Rao 
et al. 1992) . The results did not show any 
significant effect of local topography on 
within a single scene for the type of terrain in 
the study area. Multiparameter least square 
regression analysis was used to fit the soil 
moisture data with multifrequency, 
multipolarization data to identify the best 
frequency - polarization combinations for the 
accurate estimation of soil moisture. There were 4 
corn fields in the study area for which surface 
soil moisture data (0-5 cm) were available for 
July 10, 13, 15 and 17. The plant water content 
of these fields varied from 1.9 to 2.2 kg/sq.m. 
Neglecting the differences in the vegetation 
parameters, a* of the 4 corn fields was plotted 
as a function of soil moisture. It was observed 
that was linearly related to soil moisture in 
accordance with Ulaby et al. (1981). However, the 
degree of nonlinearity is more for L-band than for 
C-band. 

Table 1 gives the correlation coefficient (r) , 
standard error of estimation of soil moisture (s) 
and slope for C-, L- and P-bands; HH, W and HV 
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polarizations; both for linear and nonlinear 
fittings. It can be seen from Table 1 that C-band 
showed higher r compared to L-band, and P-band 
showed the least. The trends were similar for s 
values. L-band showed higher slope values 
compared with C-band, indicating the higher 
sensitivity of L-band to soil moisture. The 
lesser sensitivity of C-band can be understood due 
to its less penetration ability through 
vegetation. Nonlinear fitting considering a 
square term showed a marginal improvement in the 
results . 

Multifrequency linear least square fit was 
attempted considering several combinations of 
frequencies and polarizations and the results are 
shown in Table 2 . It can be seen from Table 2 
that the combination which includes all 
polarizations and frequencies gave the best 
correlation coefficient (r = 0.9) and standard 
error of estimation (s = 2.83). However, other 
combinations are also reasonably good. For 
further investigations, the combination of L-, C- 
bands with HH, W and HV polarizations was 
considered which showed r = 0.89 and s = 2.94. 

It should be noted that, in all the above 
calculations, the same data were used for 
computing the coefficients as well as for 
estimating r and s. To check if the coefficients 
are useful for other data sets, we estimated the 
coefficients using the data sets of July 10, 13 
and 15. These coefficients were then used with 
of July 17 to estimate the soil moisture. 
Comparison with field measured soil moisture 
showed good agreement between estimated and 
measured soil moisture for July 17. To assess how 
best these coefficients can be used for other corn 
fields in the study area, we estimated soil 
moisture values for some other corn fields and 
observed that they were on the same order of 
magnitude as that of the field measured soil 
moisture for all 4 dates. There were no ground 
truth data available for these fields to directly 
verify the results. However, the meteorological 
station within the study area recorded 1.3, 5.0, 
18.1, 18.5, 19.8 and 0.0 mm for the dates July 10, 
11, 12, 13, 16 and 17. The estimated soil moisture 
was in accordance with the precipitation rates. 


3. CONCLUSION 

In spite of the presence of corn fields in the 
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study area, SAR data showed good correlation with 
soil moisture for L — and C — bands. A simple linear 
regression of multipolarization, multifreguency 
approach gave correlation coefficient of 0.9 and 
standard error of estimation of 2.83. The 
regression coefficients derived from the 4 corn 
fields with ground data were useful for 
estimating soil moisture levels of other corn 
fields in the study area for the same date as well 
as for other dates. One limitation of the current 
approach is that the regression analysis 
considered each frequency meassurement as 
independent. Therefore, there is a need to 
develop a model which makes use of the inter- 
relationships among frequencies (in terms of 
penetration capability) and polarizations (in 
terms of geometry of vegetation) in such a way 
that the effects of vegetation can be accounted 
for automatically. Our ongoing research addresses 
this problem. 
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TABLE 1. RESULTS OF THE LINEAR AND NONLINEAR 
REGRESSION ANALYAIS 


DATA LINEAR REGRESSION NONLINEAR REGRESSION 

o* r s SLOPE r s 


C-HH 

0.75 

4.23 

0.17 

0.77 

4.10 

C-W 

0.83 

3.57 

0.19 

0.83 

3.56 

C-HV 

0.79 

3.96 

0.13 

0.79 

3.95 

L-HH 

0.77 

4.06 

0.32 

0.81 

3.73 

L-W 

0.61 

5.09 

0.40 

0.76 

4.16 

L-HH 

0.79 

4 . 00 

0.29 

0.79 

3 . 94 

P-HH 

0.19 

6.30 

0.16 



P-W 

0.20 

6.29 

0.15 



P-HV 

0.09 

6.69 

-0.05 




TABLE 2 . RESULTS OF LINEAR LEAST SQUARE FITTING 

FOR SEVERAL COMBINATIONS OF FREQUENCY AND 
POLARIZATION 


S.NO 

BAND 

POLARIZATION 

r 

s 

1 

C 

HH, W,HV 

0.84 

3.51 

2 

L 

19 

0.85 

3.33 

3 

P 

If 

0.31 

6.11 

4 

C,L,P 

HH 

0.80 

3.83 

5 

II 

W 

0.84 

3.44 

6 

VI 

HV 

0.84 

3.52 

7 

C, L 

HH, W,HV 

0.89 

2.94 

8 

C,P 

II 

0.86 

3.27 

9 

L,P 

II 

0.87 

3.16 

10 

C / L, P 

II 

0.90 

2.83 

11 

C , L 

HH, HV 

0.86 

3.22 

12 

C,P 

ii 

0.82 

3.68 

13 

L,P 

vi 

0.86 

3.26 

14 

C, L 

HH, W 

0.85 

3.33 

15 

C,P 

ti 

0.86 

3.30 

16 

L,P 

ii 

0.84 

3.48 

17 

C/ L, P 

HH, W 

0.87 

3.15 

18 

If 

HH, HV 

0.87 

3.17 
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This paper reports on results of ongoing efforts to develop an algorithm for 
soil moisture retrieval from SAR imagery. Estimates of soil moisture are of great 
importance in numerous environmental studies, including hydrology, meteorology, 
and agriculture. Previous studies [1] using extensive scatterometer measurements 
have established the optimum parameters for moisture retrieval as C-band HH radar 
operating at incidence angles between 10° to 15°. However, these parameters have 
not been tested or verified with imaging radar systems. The results from different 
investigators have shown considerable variability in the relationship between soil 
moisture and radar backscattering. This variability suggests that those algorithms 
are site-specific. Furthermore, the small incidence angle requirement limits the 
spatial application, especially for airborne radar systems. 

The imaging radar polarimeter permits measurement of the full polarization 
signature of every resolution element in an image. The radar polarization signa- 
ture of an object permits a more accurate description of the object of interest than 
single-polarization measurements [2]. Thus, the solution for geometric shape and 
dielectric constant of an object is less ambiguous, making the development of a quan- 
titative algorithm for soil moisture retrieval from Synthetic Aperture Radar (SAR) 
data possible. Our previous work [3] indicated that the ratio of the co-polarization 
signals, that is, the ratio of cr hh to <r vv , could be used for soil moisture retrieval at 
longer wavelengths (L-band) and at larger incidence angles ( > 40° ). The algorithm 
to infer soil moisture from imaging radar data was based on a first-order surface 
scattering model. This model predicts that the co-polarization ratio is sensitive to 
soil moisture at large incidence angles but not to surface roughness. However, the 
polarization signal ratio measurements are sensitive to the radar system noise and 
the other scattering contributions such as the multi-surface and volume scattering 
even if these effects only contribute a small portion of the total signal in the mea- 
surements. These factors result in an under-estimation of soil moisture when the 
first-order surface scattering algorithm was applied to imaging radar data. 

To address these problems, JPL AIRSAR data, in May 1988 and in September 
1989 acquired over an agricultural area near Fresno, California, were used to test 
the algorithm. We evaluated the effects of the radar system noise, the multi-surface 
scattering [4] and the volume scattering from soil [5] on the co-polarization ratio 
measurements. 

Assuming the noise power to be the same in all channels, the effect of system 
noise on the ratio measurement of <r vv /a hh can be expressed as 

a vv + noise <7 VV ^ 

cr hh + noise a hh 

for all values under the condition of cr hh < <r vv . This results in an under-estimation 
of the soil moisture especially at larger incidence angle range because the signal to 
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noise ratio typically decreases as incidence angle increases. Based on the assump- 
tion that the noise power in the two co-polarized channels is typically of the same 
magnitude, the amount of noise in all channels can be estimated. This estimated 
noise can be then used to adjust the observed measurement before inferring soil 
moisture. We will show the principle and comparison of the estimated with the 
measured system noise power. 

The surface scattering models assume that the scattering medium is a homoge- 
neous dielectric half-space. In practice, natural soil is not a perfectly homogeneous 
dielectric medium. Instead, it is a mixture of soil particles, air pockets, and liq- 
uid water. This results in dielectric discontinuities inside the soil. Because soil is a 
densely packed medium, the effects of these discontinuities will be reduced for longer 
wavelengths, especially when the distance between scatterers is much smaller than 
the wavelength. The result is that the volume scattering of soil contributes only a 
small portion of the observed signals at longer wavelengths and that the dominant 
scattering source is the surface backscattering at the air-soil interface. In evaluat- 
ing the magnitude of each co-polarization signal, the surface scattering can be used 
to explain the general relations between the backscattering measurements and soil 
physical properties. However, in attempting to relate the polarization ratio or differ- 
ence to the physical properties of soils, the volume scattering contribution becomes 
significant even if it only contributes a small portion in the observed backscattering 
returns. This effect is also expected when using long-wavelength sensors because 
of deeper penetration. To overcome the volume scattering effect on estimation of 
soil moisture, we have developed an algorithm which is based on the first-order 
scattering model considering both the surface and volume scattering contributions 

[5]- 

Figure (A) shows an image of the inferred soil moisture map of the study 
sites from L-band SAR data. This map was produced using an first-order surface 
scattering model only. The soil moisture map shown in (B) was derived by the 
algorithm which includes both surface and volume scattering of soil. The black 
regions are vegetation covered fields. When applying the first-order surface scatter- 
ing algorithm, only about 20 to 30 percent pixels were within the possible physical 
conditions predicted by the first-order surface scattering model. As shown in Fig- 
ure (A), there are many pixels with missing values even after post-processing. It is 
especially evident at large incidence angles. However, applying the algorithm with 
both surface and volume scattering considerations, about 80 percent of the pixels 
were within the physical limits. During the SAR flights, the volumetric soil mois- 
ture for the sampled dry fields varied between 3 and 10 percent. Most bare fields 
were dry because none of them had been irrigated for at least several weeks. The 
inferred soil moisture from SAR data agrees well with the field measurements and 
value ranges from 2 to 14 percent were inferred. 

To evaluate the polarimetric SAR parameters for retrieval of bare soil moisture, 
we examined (1) the applicability of the first-order surface scattering model through 
the measurements of the depolarization factor at P-, L-, and C-band from bare fields, 
and (2) all ratios of the co-polarized signals and their linear combinations through 
the model predictions. 

The measurements of the depolarization factor at C-band indicate a signif- 
icant multi-scattering involved in the C-band measurements, but not at L- and 
P-band. When the surface is relatively rough, the higher-order surface scattering 
decreases the difference between the co-polarization signals. This also causes an 
underestimation of the soil moisture if we apply the first-order inversion algo- 
rithm. When second-order terms are added in the scattering calculations [2], the 
co-polarization ratio measurements become sensitive to roughness differences, but 
this sensitivity is reduced at larger incidence angles and for smoother surfaces. 

First-order surface backscattering models predict that the ratio of the backscat- 
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tering coefficients of the co-polarizations is sensitive to soil moisture and insensitive 
to surface roughness. Similarly, the ratios of the linear combinations from the co- 
polarized returns have the same properties. The task is to select the measurements 
which maximize the sensitivity to soil moisture. Through the model simulations, 
we found that in all possible co-polarization ratios, <r jcr and <7 j{c <J ) 
provide the greatest sensitivity to soil moisture at larger incidence angles. 

As discussed above, the polarimetric SAR parameters for retrieval of bare soil 
moisture are likely to be L-band, the co-polarization ratios of a vv /<j hh and a hh /(a vv - 
<r hh ), and incidence angles above 40°. 



Figure. Comparison of inferred soil moisture map by using the first-order surface 
scattering model at top in (A), and by the algorithm which includes both surface 
and volume scattering at bottom in (B). The image brightness is proportional to 
soil moisture ranging from 2 to 14 percent by volume. 
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L Introduc tion 

The JPL Airborne S AR (AIRSAR) measures the full polarimetric properties of scene 
objects. AIRSAR data contain a complete set of magnitudes {for all linearly-polarized 
combinations [v/z., VV, HH, and CS (cross polarization)]} and phase-angle differences 
(for all pairs of linearly-polarized components) for each image pixel. Radar image users 
have been greatly interested in the like-polarization (LK) phase-angle difference 
parameter, Acp^. Usually, it is associated with the type of interaction between 

microwaves and scene-object scattering-elements. One value of A(p^ is said to be 

produced by a dominant single scattering (or reflection) event, and another value is said 

to be the result of a dominant double-scattering (or double reflection) event. These are 

often called “single bounce” and “double bounce” interaction signatures, respectively. 

There appears to be a significant amount of confusion about how the interaction 

processes produce the characteristic values of A© and about what value of A© should 

lk y lk 

be associated with each specific type. This paper attempts to address this issue. 

2, Polarized Reflection Coefficients for a Smooth Interface 

Most electromagnetic (EM) wave theory references present equations for calculating 
the (electric-field) reflection coefficients, p and p . of a smooth interface between two 

media, where p^ is for the case of the electric-field vector oscillating in the direction 

normal to the plane of incidence [viz., the TE or horizontal polarization ( H ) case], and p^ 

is for the case of the electric-field vector oscillating within the plane of incidence [v/z., 
the TM or vertical polarization (V) case]. It is important to note that p is a complex 
number which relates the magnitude and phase of the electric field of the reflected wave 
to the magnitude and phase of the electric field of the incident wave near the interface . 

There are many ways to express the p equations. One set is as follows: 

nj cos - n 2 cos 0 2 n i cos $2 " n 2 cos 

P H = and p y = (la, lb) 

nj cos 0j + n 2 cos 0 2 nj cos 0 2 + n 2 cos Q l 

where n, is the complex index of refraction for Medium 1, n 2 is the complex index of 
refraction for Medium 2, 0! is the incidence angle in Medium 1, and 0 2 is the refraction 
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angle in Medium 2. The reader may notice that Eq. (lb) has the opposite sign to 
expressions for p y given in some reference materials. In other words, the numerator 

terms in Eq. (lb) are sometimes reversed in standard references. 

It is important to understand the reason for the sign discrepancy in expressions for 
p . Interestingly, a sign discrepancy never occurs in expressions for P H > The genesis of 

the error is the way in which the reflection problem is posed. For the H case, the incident 
and reflected electric fields are both parallel to the interface and are, therefore, parallel to 
each other. Thus, phase-angle comparisons are easy to made for this case. For the V 
case, however, the incident and reflected electric fields are, in general, not parallel to each 
other (see Fig. lb). Since phase-angle comparisons between oscillating electric fields 
must be made for parallel components, one has to chose which set of components to use 
for the V case. For this case, one set of electric field components is perpendicular to the 
interface, and the other set is parallel to the interface. When the perpendicular 
components of the incident and reflected waves are chosen, the derived expression 
for p has a sign error. When the components parallel to the interface are chosen, 

the sign for p y is correct. This simple error has led to the confusion referenced above. 

3. The Single-Bounce Case 

The correctness of Eq. (lb) is easily verified by considering the case where 0i = 0 2 = 
0 (i.e., the case of normal incidence). This is the case for a single-bounce interaction that 
can produce significant “backscattering” to a SAR imaging system. In this case, both p h 

and p y are equal to [(rij - n 2 ) / (n t + n 2 )]. If the sign of Eq. (lb) were reversed, the 

expected equality of reflection coefficients for normal incidence fails to occur. Also, note 
that both reflection coefficients are always negative for the usual case where > n j. An 
important implication of the above is that the phase angles of the electric fields for 
both polarizations change by 180 degrees due to normal (single-bounce) reflection 
from a denser medium . If Medium 1 and/or Medium 2 are lossy (i.e., have significant 
imaginary parts), then the phase angle shift will differ somewhat from 180 degrees. 
Nevertheless, for normal incidence, P H *and p y are the same , and A(jy K equals 0. 

Since a double-bounce interaction appears to be simply two single-bounce 
interactions, one after the other, one might believe that A(p ) K will still be zero for this 

type of interaction. This is not the case. AIRSAR data analysts often observe values of 
Atp for double-bounce objects that are 180 degrees away for the values of Atp^ for 

single-bounce objects. The solution to this perplexing inconsistency is given in the next 
section of this paper. 

4. The Double-Bounce Case 

Consider the double-bounce case for H polarization in Fig. la. Due to the reflection 
from Interface A, the H wave undergoes a reversal of phase. For the usual case of n 2 > 
n,, (n 2 cos 0 2 ) is always greater than (n, cos 0|). Thus, the phase of the reflected H wave 


73 


will always be reversed by the reflection. Another phase reversal occurs in the H wave as 

the result of the second reflection from Interface B. Thus, the phase-angle difference 

between the incident wave and the comer-reflected “backscattering” wave will be near 

zero (depending on the material complex indices of refraction for the media on both sides 

of the two interfaces involved). The V wave (see Fig. lb) undergoes similar changes in 

phase during each of the two reflections; however, one must be careful to apply the phase 

reversals to the electric field components of the V wave which are parallel to the two 

interfaces involved. When this is done, one sees that the comer-reflected “backscattered” 

wave will be approximately 180 degrees out of phase with the incident wave. Therefore, 

the difference in the phase-angles of the H wave and the V wave, Atp , will be about 

LK 

180 degrees for this double-bounce interaction. Thus, it is the geometry of the comer 

reflection itself, that “turns” the direction of the electric field for V-polarization without a 

corresponding “turning” the H polarization, which causes the unique like-polarization 

phase-angle difference signature of double-bounce dominated scene objects [e.g., 

wetlands and some forests and woodlands (with a smooth, wet substrate)]. An exception 

will occur when any local angle of incidence exceeds the Brewster angle on either 

interface. This (Brewster condition) happens when (n, cos 0 2 ) is greater than (i^ cos 6,). 

In this case, the phase shift of the V wave is affected, and Atp reverts to near zero. For 

L.IC 

angles near the Brewster condition, the magnitude of the V component also drops towards 
zero; therefore, the use of phase-angle differences for near the Brewster condition is ill 
advised. Other propagation phenomena also affect Atp^ For example, birefringent 

media for some volume-scattering situations (e.g., a com crop with highly-oriented 
vertical stalks) will also alter Atp through the effects of differential speeds of 

propagation with polarization. Since single-scattering interactions affect more pixels than 
other interactions, one may calibrate Atp by noting the mode of its distribution. 

LK 



Figure 1A. Horizontal Polarization. 



Figure IB. Vertical Polarization. 


74 


N 94-16757 


APPLICATION OF MODIFIED VICAR/IBIS GIS TO 
ANALYSIS OF JULY 1991 FLEVOLAND AIRSAR DATA 

L. Norikane, Broek, B.*, and Freeman, A. 

Jet Propulsion Laboratory 
California Institute of Technology 
4800 Oak Grove Drive 
Pasadena. CA 91 109 

TNO-FEL 
250g JG Den Haag 
The Netherlands 


1. INTRODUCTION 

Three overflights of the Flevoland calibration / agricultural site 
were made by the JPL AIRSAR on July 3, July 12, and July 28. 1991 
as part of MAC-Europe '92, A polygon map was generated at TNO- 
FEL which overlayed the slant range projected July 3 data set. Each 
polygon was Identified by a sequence of points and a crop label. The 
polygon map was composed of 452 uniquely identified polygons 
and 15 different crop types. 

Analysis of the data was done using our modified VICAR/IBIS GIS 
[1], This GIS is an extension of the VICAR/IBIS GIS first developed 
by Bryant [2] in the 1970's which is itself an extension of the VICAR 
(Video Image Communication and Retrieval) image processing 
system also developed at JPL [3] . 

The VICAR/ IBIS GIS requires that all images be coregistered to a 
georeference image. The georeference image must contain uniquely 
identified polygons representing homogeneous areas such as 
counties or vegetation stands. After conversion to a ground range 
projection, the Flevoland polygon map was ideally suited for use as 
the required georeference image. 

To integrate the AIRSAR data into the VICAR/IBIS GIS. we 
generated 5 VICAR-labeled ground range projected Images for each 
frequency and each data set (or overflight). The image set consisted 
of HH. HV. W, HHW* phase, and HHW* amplitude images. For 
each image, we also generated an "info" file which contained a table 
of conversions from pixel value to "actual" value such as radar 
backscatter or radians. 

The images from July 12 and July 28 were coregistered to the July 3 
images using a two step process: image rotation/ scaling and image 
warping. Both steps use tiepoints which are selected by hand. 

Once all images had been coregistered, the following statistics were 
calculated for each polygon, frequency, and data set: average HH. 
HV. and W backscatter, average HHW* phase and amplitude, 
correlation coefficient using the averaged statistics, odd-bounce, 
double-bounce, and volume scattering percentage contributions [4], 
and average Incidence angle (per data set only). These data were 
stored in columns of the georeference info file. 
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For each data set, selected columns In the Info file were written Into 
a text file which served as Input to a Bayes classifier [5J. Depending 
on the features selected, the classifier identified each polygon with 
a high degree of accuracy. 

2. PRELIMINARY RESULTS 

We attempted to classify each polygon into 1 of 9 classes. (6 of the 
original 15 classes were eliminated due to insufficient numbers of 
representative polygons.) To do this we used a Bayes classifier and 
assumed Gaussian distributions. 

Different combinations of features from the July 3 and July 28 data 
sets were used to classify the polygons with varying degrees of 
success. 

Figure 1 shows the results from using only C-band HH, HV, W, 
HHW* phase, and correlation coefficients from both data sets. 

Using only these features yielded good classification results (i.e. no 
polygon was mlsclassified more than 25% of the time) 


Features used: 


July 3 / C-Band HH, HV, W. HHW* phase, correlation coefficient 
July 28 / C-Band HH, HV, W, HHW* phase, correlation coefficient 


Actual 

Percent Classified As... 






Crop 










Type 

RAP 

GRA 

POT 

WHE 

SBT 

COR 

BAR 

BEA 

FRU 

RAP 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

GRA 

0.0 

87.5 

4.2 

2.1 

4.2 

2.1 

0.0 

0.0 

0.0 

POT 

0.0 

1.0 

96.1 

0.0 

1.9 

1.0 

0.0 

0.0 

0.0 

WHE 

0.0 

4.3 

0.0 

91.5 

4.3 

0.0 

0.0 

0.0 

0.0 

SBT 

0.0 

3.8 

0.0 

0.0 

96.3 

0.0 

0.0 

0.0 

0.0 

COR 

0.0 

6.3 

6.3 

0.0 

12.5 

75.0 

0.0 

0.0 

0.0 

BAR 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

100 

0.0 

0.0 

BEA 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

100 

0.0 

FRU 

0.0 

0.0 

0.0 

0.0 

6.3 

0.0 

0.0 

0.0 

93.8 

| RAP = 

rapeseed, GRA 

= grass, POT 

= potatoes, WHE = 

wheat. 

SBT = 


sugarbeet, COR = com, BAR = barley, BEA = beans. FRU = fruit) 


FIGURE 1 


Other combinations of features have been tried. Results will be 
presented at the workshop. 

From our analysis thus far, it is clear that there are more features 
available for classification than are truly needed. Additional work 
needs to be done in order to assess which features are highly 
correlated and thus not necessary for classification. 
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3. SUMMARY 

We have successfully Integrated the three 1991 Flevoland 
calibratlon/agricultural site AIRSAR overflight data sets In our 
modified VICAR/IBIS GIS. The Images were successfully 
coregistered to a georeference Image generated from a polygon map 
which overlayed the July 3 data set. Statistics and other 
information were accumulated and saved In the georeference info 
file. Selected columns from the info file were used as features In a 
Bayes classifier assuming Gaussian distributions. Using these 
features polygons were correctly classified at least 75% of the time 
In the example shown. Additional work needs to be done to assess 
which features are highly correlated and thus not necessary for 
classification. 
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1. INTRODUCTION 

Imaging radar data provides information about the geometric and dielectric 
properties of the Earth's surface (Ulaby et al, 1982). The Jet Propulsion Laboratory (JPL) 
polarimetric Airborne Synthetic Aperture Radar (AIRSAR) is currently obtaining imaging 
radar data for use in geologic, vegetation, snow and ice, and oceanic studies (Held et al., 

1988; Zebker et al, 1987). In the near future, the Shuttle Imaging Radar C (SIR-C/X- 
SAR) and the Earth Observing System Synthetic Aperture Radar (EOS SAR) will also 
collect polarimetric radar data (JPL, 1 986, 1991). 

A need exists for a user-friendly, interactive software package for analysis of 
these polarimetric radar data sets. Previous software packages, such as JPL's Multiview 
(Burnette, 1988), while providing some analysis capabilities for these data, did not allow 
interactive viewing and were tied to specific image display hardware with operating 
system dependencies. 

A prototype software system, the "Radar Analysis and Visualization 
Environment" (RAVEN) under development at the Center for the Study of Earth from 
Space (CSES) at the University of Colorado, is designed to deal with data from the JPL 
AIRSAR instrument and other proposed polarimetric radar instruments. RAVEN is being 
developed using the Interactive Data Language (IDL) (Research Systems, 1991). It takes 
advantage of high speed disk access and fast processors running under the UNIX operating 
system in a X-windows environment to allow for rapid, interactive visualization of 
AIRSAR data in both image and graphical ways. It provides a user-friendly interface 
through the use of menus, sliders, buttons, and display windows (Figure 1). 

2. DISPLAY FEATURES 

Analysis of polarimetric radar data requires innovative display of both images and 
quantitative information. The images provide spatial information and the graphics 
provide frequency and polarimetric information. It is necessary to view an entire scene 
and at the same lime be able to see subscencs at full resolution, and enlarged to observe 
details. RAVEN is designed to create and manage multiple data windows. An entire 
AIRSAR scene subsampled to fit in a scroller window allows selection of an area to be 
displayed at full resolution in a image window. A zoom window allows for a subset of 
the image to be panned and zoomed from 1 to 16 times. RAVEN also provides the 
capability to interactively select and display synthesized polarimetric images for three 
frequencies as gray scale images, density-sliced images, or as color composites. These 
images are displayed as sigma-zero backscatter or total power values. Images of any 
polarization combination can be synthesized in RAVEN using the Stokes matrices for the 
calculations. This synthesis, however, takes approximately 2 minutes and is not 
considered real-time. Therefore, RAVEN also uses a band sequential format data set that 
contains various pre-synthesized polarization images. This makes possible rapid display 
of gray scale and color images. 
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Rapid viewing of frequency and polarization information for individual pixels 
allows users to get an idea of the variation in a data set. It also allows definition of areas 
for detailed analysis. RAVEN uses co-polarized and cross-polarized signature windows 
that provide real-time viewing of cross sections of polarization signatures. A window is 
also provided to view frequency plots. Frequency plots and full polarization signatures 
from multiple pixels are saved in another window and allow the user to compare different 
units. 

3. ANALYSIS FEATURES 

The analysis capabilities of RAVEN presently include creating and viewing 
polarization and frequency signatures for the three AIRS AR frequencies using a pull- 
down menu in the main RAVEN window. Calculation of polarization signatures and 
frequency plots for both individual pixels and polygon regions is provided. Polarization 
signatures are 3-dimensional plots that represent the radar backscatter as a function of all 
ellipticity and orientation polarizations (van Zyl et al, 1987). The polarization signatures 
generated in RAVEN are displayed as both 3-dimensional plots and 2-dimensional images. 
These plots and images are displayed as either gray scale or density-sliced color. 
Polarization signatures extracted from the data can be compared to each other or to 
theoretical polarization signatures created for specified input scattering matrix values. A 
signature viewer window also allows both subtraction and division operations to be 
performed on polarization signatures. 

Images of RMS surface roughness, fractal dimension, and dielectric constants are 
imported from an inversion program and displayed in RAVEN (Kierein-Young and Kruse, 
1992). The images allow the user to view the spatial variance of these parameters and 
help in creating maps of different units. 

4. FUTURE DEVELOPMENT 

RAVEN is a prototype software tool under development for analysis of 
polarimetric radar data. Features are being added to deal specifically with analysis 
requirements for geologic targets. Classifications, library functions, polarimetric 
matching techniques, and inversion models are some of the future enhancements planned 
for the RAVEN package. 

5. ACKNOWLEDGEMENTS 

Portions of this work were supported under a NASA Graduate Student Researchers Grant 
(NGT-50728). Portions of this work were supported under the SIR-C/X-SAR project by 
NASA/JPL contract 958456. 

6. REFERENCES 

Burnette, Fred, 1988, JPL Multipolarization Software Package, JPL internal publication, 
142 p. 

Held, D. N. et al. 1988, The NASA/JPL multifrequency, multipolarization airborne SAR 
system: in Proceedings. IGARSS '88 . The Institute of Electrical and Electronics 
Engineers, Inc. (IEEE), New York, p. 345-349. 

Jet Propulsion Laboratory, 1986, Shuttle imaging radar-C science plan: JPL Publication 
86-29 . Jet Propulsion Laboratory, Pasadena, California. 

Jet Propulsion Laboratory, 1991 , The Earth Observing System Synthetic Aperture Radar 
(EOS SAR): Report to the Space Science and Applications Advisory 


79 


Committee (SSAC) for the Office of Space Science and Applications (OSSA) 
1992 Strategic Plan, 25 p. 

Kierein- Young K. S., and Kruse, F. A., 1992, Extraction of quantitative surface 

characteristics from AIRSAR data for Death Valley, California: in Proceeding s 
of the Fourth AIRSAR Workshop (this volume). 


Research Systems, Inc. (199 1Y IDL® User's Guide Version 2.2 . 

Ulaby, Fawwaz T„ Richard K. Moore, and Adrian K. Fung, 1982, Microwave Remote 
Sensing Active and Passive . Addison-Wesley Publishing Company, 
Massachusetts. 


van Zyl, J. J., H. A. Zebker, and C. Elachi, 1987, Imaging radar polarization signatures: 
Theory and observation, Radio Science . V. 22, No. 4, pp. 529-543. 

Zebker, Howard A., Jakob J. van Zyl, and Daniel N. Held, 1987, Imaging Radar 

Polarimetry From Wave Synthesis, Journal of Geophysical Research . V. 92, 
No. Bl,pp. 683-701. 



I . . - ' • ' • ■ 

a T * S * V ‘ Pfpiu y Sl1 ^ f Uf 1 ny Scale [ ... 


Figure 1 . An example of the RAVEN software display screen. 
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ABSTRACT 

Using the JPL AIRSAR interferometer, measurements of the ocean coherence time at L and C band 
can be made at high spatial resolution. Fundamental to this measurement is the ability to image the 
ocean interferometrically at two different time-lags, or baselines. By modifying the operating 
procedure of the existing two antenna interferometer, a technique has been developed make these 
measurements. L band coherence times are measured and presented in this paper. 


Introduction 


Over the past 20 years or so there have been many theories and models put forth concerning 
microwave sensing of the ocean surface and the importance of the time scale at which the scattering 
remains coherent. The coherence time in SAR is particularly important since SAR imaging relies 
on coherent integration to achieve the along-track resolution. Some recent estimates of the ocean 
coherence time at L band have been made on the order of 0.05 to 0.1 seconds [1,2,3J. 


The along-track interferometric SAR used by JPL [4] to measure ocean surface velocities makes 
a direct measurement of the coherence of the ocean surface over the lag time separating the 
observations by the two antennas. A technique has been devised [5] where measurements at two 
baselines, or lag times, may be acquired simultaneously with the existing AIRSAR two antenna L 
band system by taking advantage of the multichannel transmit and receive capability required tor 
polarimetric data acquisition. The AIRSAR polarimeter acquires full polarization information by 
transmitting H and V alternately while receiving H and V from each transmit event. Using an 
identical scheme, radar pulses are emitted from the forward (F) and aft (A) antennas alternate y 
and received by both. This produces four channels of data, FF, FA, AF and AA, which may be 
processed to imagery analogous to the polarimetric case. Two unique interferometric baselines 
may be constructed using these data, one being the separation distance of the antennas and the 
second being half this distance. 


Example data 


The AIRSAR instrument imaged the Strait of Messina in the summer of 1991. The mode of 
operations for one data take was the dual ATI baseline mode for L and C, and standard quad-pol 
for P band. A phase difference image is obtained (after processing to remove aircraft motion and 
sensor geometry phase effects) by combining the LAA complex image (A) with (in this case) the 
LAF (F) complex image according to the following formula: 

A<p = tarr 1 ( l ™eal(C ^ (1) 


where: 

C - F A* 


( 2 ) 
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A multi-look phase is determined by coherently summing the interferogram pixels, C, and 
retaining the resultant phase. This method has been shown by Rodriguez [6] to be the maximum 
likelihood estimator for determining the interferometric phase. In this case summation over 16 
azimuth pixels was used. The phase difference may be converted to interferometric velocities 
using the following equation: 


u = A<p 


Av 

4nB 


(3) 


where A(p is the observed phase change, B is the baseline length, v is aircraft speed and A is the 

radar wavelength. The ratio of B/v represents the time lag, T, between forward and aft imaging. 
The interferometer velocity represents the radial component of the phase change. To project this 
onto the ground image plane, to extract the surface component, one must divide by sin(Q), where Q 
is the incidence angle. The surface interferometer velocity may be written as: 


u s = 


A(p 
kb * 


(4) 


where kb = 4nlX sin( 6) and is identified as the Bragg wavenumber. 

Measuring coherence time with ATI 

In addition to observing the phase difference and producing interferometer velocity maps, 
correlation maps of the ocean surface may be produced for each baseline. With this dual baseline 
information, the coherence time of the ocean surface can be measured directly. The coherence time 
is defined here as the time it takes the autocorrelation function of the backscatter field from the 
ocean to be reduced to He. The correlation strength between two images forming a baseline is 

given by p,- where i denotes which of the two baselines the correlation calculation is being carried 
out on. The correlation image may be calculated as follows: 

Z/a * 

V 11/ 1 2 51 a I 2 

where/ and a are the pixels of the forward and aft phase center images for baseline i, and the 
summations are carried out over an adjacent number of pixels, typically 8 or 16, similar to a 
multilook filter. 




By using the two baselines available with the AIRSAR data, one may measure two values of p: p s 

and pi where s and / indicate “short” and “long” baselines. The nominal time lags associated with 
the two AIRSAR baselines corresponding to these measurements are 0.05 and 0.10 seconds 

respectively. Assuming a Gausian correlation function, the lie coherence time, t c> may be 
calculated using these two measurements: 





Ts 2 + ?/ 2 
In (psPl) 


( 6 ) 
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High resolution correlation images for the short and long baseline of the Messina data have been 
formed using (5). From these short and long baseline correlation images, a coherence time image 
has been formed using (6). The image was quantized to 16 levels going between 0.0 seconds and 
0.32 seconds in linear units of 0.02 seconds. The land areas exhibit coherence times greater than 
the 0.32 maximum. There is quite a bit of structure to this map, both large scale and small scale. 
There are areas in which the coherence time is noticeable larger than other areas. Ocean coherence 
times observed in this L-band image vary from about 0.09 seconds to 0.27 seconds. 
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In response to mounting scientific evidence that the biosphere and atmosphere of 
the earth are undergoing significant anthropogenic changes that appear to be increasing in 
magnitude and rate, the United States launched the Global Change Research Program 
(USGCRP) in 1989. The goal of the program is to "...Gain a predictive understanding of 
the interactive, physical, geological, chemical, biological and social processes that 
regulate the total earth system.” Seven government agencies are involved in the program 
with the National Aeronautics and Space Administration (NASA) playing a key role by 
providing space based and earth based data. 

The mission elements of NASA’s approach to global change research involve the 
collection and analysis of satellite data from earth probes, geostationary platforms and the 
primary element of the NASA effort, the Earth Observing System (EOS). EOS consists 
of a series of polar orbiting platforms that will carry a variety of sensors capable of 
collecting digital data with spectral and spatial resolution suitable to support multi- 
disciplinary research in earth system science. 

Many of the new remote sensing systems that will begin generating data by the 
end of this century are based on experiments conducted, and knowledge gained, using 
airborne platforms, especially within the last 25 years. Although the capability to conduct 
large scale monitoring of earth processes may be dated from July, 1972 with the launch of 
Landsat 1, neither the Landsat system nor any of the numerous devices launched since 
Landsat 1, or now planned for launch, preclude the need to develop and operate remote 
sensing systems on aircraft. Airborne remote sensing will continue to have scientific 
merit because 1) new concepts in remote sensing can be developed and proved faster and 
more cheaply on aircraft than on satellites; 2) the capabilities of aircraft systems are likely 
to improve, e.g. long duration, pilotless flights, increasing the utility of aircraft for earth 
observation measurements, especially observations required in response to unusual, 
rapidly developing events; 3) measurements from aircraft have been shown to be critical in 
studies of the expression of biophysical processes in the landscape; and 4) evidence 
indicates that data from airborne sensors increases the value of satellite derived data on 
global change processes by supplying information that lowers the error associated with 
satellite-data-driven estimates of the output of those processes. 

A bibliography of global change, airborne science from 1985 to 1991 was 
assembled that may be of utility to the science community. The bibliography was 
compiled using EndNote® Plus, a software system that operates in Macintosh and PC 
environments. The intent of the bibliography is to include all articles (primarily in 
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refereed journals) that address global change issues and that utilize aircraft to acquire data. 
The bibliography is a tool that has a variety of functions. It can be used by the science 
community as a reference to search for articles on specific topics, and it is gauge for 
assessing the impact of airborne sensors on global change science. From the bibliography 
data base information can be extracted on who is performing global change airborne 
science, who is sponsoring the research, where the research is being performed and what 
type of research is being done. 

The bibliography includes entries for journal articles, conference proceedings, book and 
book sections, and reports. For journal articles the data base includes author(s), title, 
journal name, volume number, year of publication, month of issue, affiliation of 
author(s), key words, a label field, and abstract (if available). The software allows for 
inclusive and exclusive searches of the bibliography by character string. Citations can 
formatted in numerous ways for entry as in-text citations or bibliographic references. The 
current version of the bibliography will be updated periodically, and new versions should 
be available bi-annually. 
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Soil conservation programs are becoming more important as 
the growing human population exerts greater pressure on this non- 
renewable resource. Indeed, soil degradation affects approximately 
10 percent of Canada’s agricultural land (Dumanski et al., 1986) 
with an estimated loss of 6,000 hectares of topsoil annually from 
Ontario farmland alone (Agricultural Institute of Canada, 1980). 
Soil loss not only affects agricultural productivity but also decreases 
water quality and can lead to siltation problems. Thus, there is a 
growing demand for soil conservation programs and a need to 
develop an effective monitoring system. 


Local topography, soil type, plant residue, vegetative cover, 
and surface roughness interact to govern soil erosion by wind and 
water. In general, smooth bare sandy soils are most susceptible to 
erosion. Contour farming, increasing crop residue and deep tillage 
to create enough surface roughness to break wind/water flow 
velocities are management techniques being used to reduce erosion 
risk. Topography and soil type information can easily be handled 
within a geographic information system (GIS). Information about 
vegetative cover type and surface roughness, which both experience 
considerable temporal change, can be obtained from remote sensing 
techniques. Thus all the components for a monitoring system are 
available. 


Previous research has shown that useful information for soil 
conservation applications can be extracted using visible and near 
infra-red wavelengths (Gausman et al., 1975). Atmospheric effects 
and winter snow cover over much of the northern United States and 
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southern Canada limit the data availability in this region of the 
electromagnetic spectrum. Microwave systems, both active and 
passive, provide a reliable alternative. Recent research has 
demonstrated that information about residue/plant cover and 
surface roughness can be obtained using these types of sensors 
(Ulaby et al., 1986). For example, Brisco et al. (1991) found that C- 
Band like polarization was effective for discriminating between 
different tillage classes on wheat stubble fields as a function of 
surface roughness. 


For further development of the technology to produce an 
operational soil conservation monitoring system an experiment was 
conducted in Oxford County, Ontario which investigated the 
separability of fall surface cover type using C-Band SAR data. This 
experiment obtained quantitative information on surface roughness 
and crop residue for about 30 fields in the study site and qualitative 
inf ormation on about another 250 fields. 


To date only C-HH narrow swath data (6 meter resolution, 
45-76 degrees incidence, 7-look amplitude format) have been 
analyzed. Transformed divergence separability statistics show that 
corn stubble can be separated into two classes (high and low 
stubble) and bare smooth fields can be differentiated from deep 
tilled fields. Confusion between bare fields and pasture or hay 
fields, which are also dark on the SAR imagery, could be removed 
with complementary optical data acquired any time earlier in the 
growing season. Alternatively, a GIS cover type layer could be 
generated from historic data as pasture and hay fields are recurring 
cover types on a yearly basis. This approach is proving useful for 
other agricultural classification purposes (Bedard et al., 1992). Corn 
fields that have not been harvested are also highly separable, as 
previous work has demonstrated (Brisco and Protz, 1980). 


Transformed divergence statistics will be calculated for the 
other polarizations (W, HV) and evaluated for separability of these 
conservation classes. Furthermore the analysis will be expanded to 
include multi-polarization combinations, co-polarization ratios, and 
cross-polarization ratios to see if improvements in class separability 
of conservation farming classes can be obtained with multi- 
polarization C-band SAR data. The results of these analyses will be 
related to upcoming spacebome SAR sensors to evaluate the 
potential for this application using data which should be available in 
the next couple of years. GIS models incorporating local topography 
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and soil type to create an erosion potential map are also being 
developed. Synergism of the SAR data with visible and near-infrared 
channels from SPOT and TM will also be evaluated. The 
presentation at the AIRSAR workshop will summarize the results 
that have been obtained by the time of the workshop. 


This research will continue with experiments already planned 
using ERS-1, SIR-C, and Radarsat data for further development of 
this application. It is envisioned that operational soil conservation 
monitoring systems are possible by Radarsat, or at the least by the 
Earth Observing System (EOS), timeframe. 
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1. INTRODUCTION 

Edge detection would appear to be a crucial tool for analysing 
multi-polarised, multi-frequency and multi-temporal Synthetic Aperture Radar 
(SAR) images. Edge structure provides a simple means for comparing different 
polarisations and frequencies, and for detecting changes over time. Due to the fact 
that edges and segments (homogeneous regions) are dual concepts, edge detection 
has an important role to play in identifying segments within which mean 
backscatter measurements for use in image classification can be made ( White 

1991). 

As part of a general investigation into edge detection in SAR imagery, an 
initial investigation has been carried out into the detectability and nature of edges 
in multi-polarised and multi-frequency SAR images. The contrast ratio (CR) 
operator was used to detect edges. This operator has previously been shown to 
perform well at detecting edges in single-polarised and single-frequency SAR 
images (Touzi et a.1. 1988; Oaves ei al. 1992). 

2. THE CONTRAST RATIO OPERATOR 

At each point in an image the contrast ratio across a n x n window 
centred at that point is calculated by splitting the window into two equal sections, 
estimating the mean value of the pixels in each section and taking the ratio of 
these two mean values. Ratios are calculated for splits along the horizontal, 
vertical and two diagonal orientations. The largest ratio gives a measure of the 
contrast across the window. The larger the ratio, the greater the contrast, and the 
greater the likelihood that an edge exists at that point. 

The contrast ratio image is thresholded to produce an edge map. So long 
as speckle is a stationary multiplicative noise process the false alarm rate in the 
edge map will be constant. Additionally, if it is assumed that speckle in an 
intensity image is uncorrelated and gamma distributed, the distribution of the 
contrast ratio can be determined and used to set a threshold for a desired 
probability of false alarm. The effects on the false alarm rate of known system 
induced correlation can be removed by sub-sampling within the processing 
window. 

3. APPLICATION TO MAESTRO IMAGES 

The CR operator was applied to C and P band MAESTRO single-look 
intensity images of Felt well, UK. For each band the binary edge maps produced 
from two different polarisations were compared by overlaying them; the symbol 
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“o” is used to represent this binary operation (i.e. HV o VH = HV x 2 + VH). 
This was done for ail possible combinations of two different polarisations. Pixel 
values in the resulting 2-bit images indicated: 0 - no edge at either polarisation, 1 
or 2 - an edge at one polarisation but not the other and 3 - an edge at both 
polarisations. The process was repeated with edge maps produced using different 
window sizes n = 3, 5, 17 and probabilities of false alarm p = 10~ 3 , 10" 4 , 10 -6 
(assuming exponential speckle). The correlation length in both range and azimuth 
was found to be two pixels; 2x2 sub-sampling was used within processing 
windows to remove its effect. 

4. RESULTS 

In each set of edge maps produced from the C band image using a 
particular window size and probability of false alarm, the greatest degree of 
similarity occurred between the edge maps of the cross-polarised terms (most 
values in HV o VH were 0 or 3). A lesser degree of similarity occurred between 
the edge maps of the co-polarised terms. The most significant differences occurred 
between the edge maps of a co and a cross-polarised term. However due to the 
similarity of the edge maps of the cross-polarised terms, the differences occurring 
between the edge maps of a co-polarised term and the HV term were similar to 
those occurring between the edge maps of the same co-polarised term and the VH 
term. 


These traits can be observed in figures 1 and 2 which show the results of 
edge detection using a 9 x 9 window with a probability of false alarm of p = 10~ 4 . 
Figure 1 shows the edge maps of the HH, HV and VV terms and the result of 
combining the edge maps of all four polarisations. The edge maps of the single 
terms are all different from each other. Differences arise not just from edges being 
better defined in some polarisation than in others but also from features present in 
some polarisations being completely absent in others. The edge map of the HV 
term contains significantly more detail than that of the HH term which in turn 
contains slightly more detail more than that of the VV term. The combined edge 
map is significantly more detailed than that of any single polarisation; each 
polarisation contributes something towards it. 

Figure 2 shows histograms of combinations of edge maps of the (a) HH 
and VV, (b) HV and VH, (c) HH and HV, and (d) VV and HV terms using 
the “o” operator described above. The HH o VV histogram shows that the 
percentage of edge pixels detected in the HH term (i.e. those pixels with values 2 
and 3) was greater than the percentage detected in the VV term (pixels with 
values 1 and 3). The HV o VH histogram shows that the small proportion of 
differences between the HV and VH terms is equally split between edges detected 
in the HV term and not in the VH term and edges detected in the VH term and 
not in the HV term. The HH o HV and VV o HV histograms show that the 
differences between co and cross-polarised terms is primarily due to a greater 
percentage of edges being detected in cross-polarised terms than in co-polarised 
terms. The total percentage of edges pixels in a given combination is indicated by 
subtracting the percentage of pixels with value 0 from 100. The greatest number 
of edges is detected by combining a co and a cross-polarised term and the least 
when the edge maps of the co-polarised terms are combined. 

Apart from the HV term, the P band image was much sharper than the 
C band image. This resulted in a much larger percentage of edges being detected. 
Once again the edge maps of cross-polarised terms were similar whilst edge maps 
of co-polarised terms were significantly different from those of cross-polarised 
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terms. 


5. CONCLUSION 

The CR operator has detected edge information useful for comparing 
differently polarised terms and for use in further processing. Significant differences 
exist between the edges detected in different polarisations. The most significant 
differences occurred between co and cross-polarised terms. The differences were 
small between cross-polarised terms, as reciprocity would predict. The images 
used were uncalibrated. However, assuming a linear system model, it can be 
shown that, to a good approximation calibration of this data involves only linear 
scalings of the intensity images (Quegan 1992). Such scalings would not affect the 
outcome of applying the CR operator. The fact that the images were not of 
comparable quality places more serious constraints on interpreting the results in 
terms of fundamental scattering properties of edges. 
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Figure 1. Edge maps of a MAESTRO C band image of 
an agricultural area produced using the contrast ratio 
operator with a 9 x 9 window and a probability of 
false alarm of 0.0001: (a) HH. (b) HV, (c) VV and 
(d) HH + HV + VH + HH. 
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1 . INTRODUCTION 

Radar backscatter intensity as measured by calibrated synthetic aperture radar 
(SAR) systems is primarily controlled by three factors: local incidence angle, wavelength- 
scale roughness, and dielectric permittivity of surface materials. Radar observations may 
be of limited use for geological investigations of surface composition, unless the 
relationships between lithology and the above characteristics can be adequately understood. 
In arid terrains, such as the Southwest U.S., weathering signatures (e.g. soil 
development, fracturing, debris grain size and shape, hill slope characteristics) are 
controlled to some extent by lithologic characteristics of the parent bedrock. These 
textiual features of outcrops and their associated debris will affect radar backscatter to 
varying degrees, and the multiple-wavelength capability of the JPL AIRSAR system 
allows sampling of textures at three distinct scales. 

Diurnal temperature excursions of geologic surfaces are controlled primarily by 
the thermal inertia of surface materials, which is a measure of the resistance of a material 
to a change in temperature (Gillespie and Kahle, 1977). Other influences Include albedo, 
surface slopes affecting insolation, local meteorological conditions and surface emissivity 
at the relevant thermal wavelengths. To first order, thermal inertia variations on arid 
terrain surfaces result from grain size distribution and porosity differences, at scales 
ranging from micrometers to tens of meters. Diurnal thermal emission observations, 
such as those made by the JPL TIMS airborne instrument, are thus influenced by 
geometric surface characteristics at scales comparable to those controlling radar 
backscatter. 

This paper is a preliminary report on a project involving a combination of field, 
laboratory and remote sensing observations of weathered felsic-to basaltic volcanic rock 
units exposed in the southern part of the Lunar Crater Volcanic Field, in the Pancake 
Range of central Nevada (Snyder el al, 1972). We focus here on the relationship of radar 
backscatter cross sections at multiple wavelengths, apparent diurnal temperature 
excursions identified in multi-temporal TIMS images, surface geometries related to 
weathering style, and parent bedrock lithology. 

2. GEOLOGIC SETTING 

The Lunar Crater portion of the Pancake Range consists of normal-faulted blocks 
of Tertiary rhyolitic to andesitic volcanic rocks, overlain in places by a Quaternary 
alkaline basalt field, which includes lava flows, cinder cones and maar craters (Scott and 
Trask, 1971; Snyder et al, 1972). The area described in this paper includes the southern 
portion of the basalt field, consisting of two lava flow complexes and two associated 
cinder cones. The Tertiary volcanic rocks in the area include andesitic and quartz lalitic 
lavas, and a number of distinct welded and bedded tuff units, displaying varying 
phenocryst and lithic contents and devitrification textures (Snyder et al, 1972). 

Basalt outcrops in the area show relatively mature upper weathered surfaces, with 
well-developed desert pavements consisting of basalt cobbles 1-10 cm in diameter resting 
on an aeolian-derived sediment of up to several m thickness. Felsic outcrop morphologies 
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arc highly variable, ranging from low rounded hills with relatively thick debris mantles to 
cap-rocked mesas with fresh bedrock exposures and over-steepened cliff faces as high as 
250 m. Alluvial deposits are predominantly sandy gravels, although several washes 
contain clasts up to large boulder size, reflecting high-energy depositional environments. 
Relationships at basalt flow margins indicate that these high valleys have been 
predominantly in a downwasting regime since the basalt eruptions. 

3. REMOTE SENSING DATA 

AIRSAR data for the southern Lunar Crater area includes three passes along the 
same azimuth, providing three incidence angles for many of the surfaces of interest. 

TIMS coverage consists of three sets of overlapping paired passes, acquired at 
approximately 4 a.m., 8 a.m. and 2 p.m. local time. Figure 1 shows a segment of the 
afternoon and pre-dawn Band 5 (10.3- 1 1 . 1 pm) TIMS passes, along with an L-band HH 
polarization AIRSAR image of approximately the same area. 

Changes in contrast relationships among surfaces between the afternoon and pre- 
dawn TIMS images are almost certainly related to thermal inertia variations. Comparison 
with the 8 a.m. image taken on the same day (not shown) indicates that differential 
heating related to slope/insolation effects is minimal by the late afternoon, and is 
certainly of negligible importance for the pre-dawn data. An interesting reversal in 
contrast is seen between adjacent outcrop areas of welded tuff (indicated by numbers 1 and 
2). The southern outcrop (1) shows relatively little difference in emittance from the 
surrounding alluvial surfaces in both scenes, while the northern outcrop (2) is much 
darker in the afternoon scene and much brighter in the pre-dawn scene. Although these 
data have not yet been calibrated to extract temperature or thermal inertia estimates, 
qualitative examination indicates that the northern outcrop has a higher thermal inertia 
than the southern outcrop, expressed as a smaller excursion in temperature relative to the 
surroundings. AIRSAR data for these two surfaces suggest a relationship between 
outcrop morphology, surface roughness and thermal inertia. Figure 2 shows the behavior 
of radar cross section (sigma zero) as a function of wave number (2nfk) for several 
surfaces. The northern tuff outcrop exhibits high sigma zero values at all three AIRSAR 
wavelengths, while the southern outcrop has much lower values, and a strong dependence 
on wavelength. Field observations indicate that the northern outcrop is dominated by 
fresh exposures and high local relief, while the southern outcrop is pervasively fractured at 
the cm scale, leading to low relief and a mantle of cm-scale gravel debris. The diurnal 
temperature contrast is clearly related to the presence of the high thermal inertia rock 
exposures in the north and the low thermal inertia gravel debris in the south. Similarly, 
the AIRSAR signatures reflect the steep topography in the north, which leads to high 
sigma zero values independent of wavelength. On the southern outcrop, the surface 
appears smooth at P band and relatively rough at C-band, consistent with field 
observations of the roughness scale of the gravel debris. Although mapped by Snyder et 
al (1972) as a single unit, the two outcrops display dramatically different radar and thermal 
signatures related to weathering style. Further sample analysis should illuminate the 
lithological basis for the contrasting morphology. 

At the "boulder field" site (3 in Fig.l), the presence of m-scale and larger 
boulders is probably responsible for the apparently higher thermal inertia (note bright 
"warm" signature in pre-dawn scene) of this deposit compared to the surrounding 
alluvium. This distinct roughness signature is reflected in the S AR data (bright patch in 
Fig.l), particularly at the longer wavelengths. 

The "friable tuff' exposure (4 in Fig.l) consists of rounded "towers" of erodible 
fine-grained, possibly volcaniclastic tuff. Extensive patches of sandy debris shed from 
these exposures give the area an overall low thermal inertia, and an extremely smooth 
signature at C and L AIRSAR wavelengths. The unique radar and thermal signature of 
this unit is directly related to grain size and cementation, properties related to lithological 
and (volcanic or sedimentary) depositional parameters. 
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Figure 1. a) TIMS band 5 (10.3-1 1.1 mm) image, afternoon and b) pre-dawn, c) AIRSAR L-band HH, illumination from the 
left. Scene dimensions are approximately 3x5 km. Locations marked in b are 1) south tuff outcrop, 2) north tuff outcrop 
3) boulder field, 4) friable tuff. 



4. 


CONCLUSIONS 


Preliminary analysis of radar, thermal and field observations of the southern 
Lunar Crater region indicates that textural characteristics that dominate radar and diurnal 
thermal signatures of surfaces are related (though probably not uniquely related) to the 
primary outcrop lithology. Application of AIRSAR data may be of use in analysis of 
TIMS data, in which grain size and compositional influences on emittance must be 
separated. Such combined multi-sensor and field studies also provide valuable experience 
for analysis of purely "remote" datasets, such as Magellan and Mars Observer. 
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Figure 2. AIRSAR HH polarization backscatter cross section (sigma zero) as a function 
of wave number (2 n/X). Corresponding wavelengths are: P-band 68 cm, L-band 24 cm, 
C-band 5.6 cm. Flat spectrum of north tuff indicates topographic control. Steep 
spectrum of south tuff indicates low "power" at larger scales of surface roughness. 
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